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ABSTRACT 

The  limiting  factor  affecting  the  performance  of  most  airborne  radar  sys¬ 
tems  in  detecting  targets  on  or  near  wind-swept  surface  of  the  sea  is  return 
echoes  from  sea  surface  known  as  sea  clutter.  This  report  addresses  the  clutter 
returns  from  a  maritime  environment  for  airborne  L-band  VV-polarised  radar. 
Modelling  of  radar  sea  return  is  discussed  with  consideration  to  propagation 
under  standard  atmospheric  conditions.  A  composite  sea  clutter  model  has 
been  developed,  which  is  based  upon  the  concept  of  a  two-scale  sea  surface 
model  with  directional  sea  spectrum  and  simple  facet  specular  return  at  near¬ 
normal  incidence.  Several  sea  clutter  models  have  been  collated  with  data 
collected  by  the  Naval  Research  Laboratory  four- frequency  radar  system  un¬ 
der  varying  sea  conditions,  and  polarimetric  synthetic  aperture  radar  images 
of  the  North-West  of  Australia.  The  results  of  the  comparison  have  shown  that 
the  composite  sea  clutter  model  better  predicts  the  backscattering  coefficient, 
cr°,  of  sea  surface  returns  for  varying  sea  conditions,  radar  look  directions,  and 
grazing  angles.  In  addition,  a  radar  backscattering  model  for  bare  soil  showed 
good  agreement  between  modelling  results  and  POLSAR  data  for  returns  from 
coastal  saline  flat  surfaces  surrounding  the  North-West  coastal  areas. 
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Modelling  Airborne  L-band  Radar  Sea  and  Coastal  Land 

Clutter 


EXECUTIVE  SUMMARY 

This  report  represents  the  result  of  a  study  performed  under  the  task  AIR99/047- 
“AEW&C  SUPPORT”.  The  performance  of  an  L-band  Airborne  Early  Warning  and  Con¬ 
trol  (AEW&C)  radar  system  in  detecting  targets  in  a  maritime  environment  is  significantly 
affected  by  clutter  backscatter  from  the  sea  surface  and  coastal  land  areas.  A  thorough 
understanding  of  these  return  signals  is  vital  in  assessing  the  performance  of  the  radar  in 
detecting  targets  over  sea  and  site-specific  coastal  land  areas.  An  essential  ingredient  in 
this  development  is  the  need  to  accurately  model  sea  surface  and  site-specific  coastal  land 
areas  reflection  or  backscatter  signals. 

The  primary  objective  of  this  study  is  to  develop  an  L-band  VV-polarised  sea  clutter 
model  to  quantitatively  predict  the  radar  sea  backscattering  coefficient,  a°,  as  a  function 
of  radar  parameters,  and  sea  surface  conditions  under  standard  atmospheric  conditions. 
This  report  looks  at  several  semi-empirical  sea  clutter  models  and  their  relation  to  mod¬ 
elling  backscattering  from  the  sea  surface  for  varying  sea  conditions,  grazing  angles  and 
radar  look  directions  for  an  L-band  radar.  In  addition,  a  second-order  composite  sea 
clutter  model  has  been  developed,  which  is  based  upon  the  concept  of  the  two-scale  sea 
surface  model  to  predict  the  L-band  radar  sea  backscattering  coefficient,  a°,  for  varying 
sea  conditions,  grazing  angles  and  radar  look  directions. 

These  semi-empirical  sea  clutter  models  and  the  second-order  composite  sea  clutter 
model  have  been  collated  with  data  collected  by  the  Naval  Research  Laboratory  four- 
frequency  radar  system  under  varying  sea  conditions  and  polarimetric  synthetic  aperture 
radar  images  of  the  North-West  of  Australia.  The  results  of  the  comparison  have  shown 
that  the  second-order  composite  sea  clutter  model  better  predicts  the  backscattering  coef¬ 
ficient,  <7°,  of  sea  surface  returns  for  a  wide  range  of  grazing  angles,  radar  look  directions 
and  sea  conditions.  The  composite  sea  clutter  model  also  showed  good  agreement  between 
the  modelling  result  and  POLSAR  data  for  shallow  water  returns  from  the  coastal  water 
surface  surrounding  the  Darwin  area  for  grazing  angles  between  25°  to  65°. 

The  secondary  objective  of  this  study  is  to  develop  an  L-band  VV-polarised  land  clutter 
model  to  quantitatively  predict  the  radar  backscattering  coefficient,  a° ,  from  coastal  land 
areas.  The  radar  backscattering  model  developed  for  predicting  microwave  scattering 
from  bare  soil  surface  showed  good  agreement  between  its  modelling  results  and  measured 
backscatter  values  from  coastal  saline  flats  surrounding  the  North-West  coastal  areas. 
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1  Introduction 

The  limiting  factor  affecting  the  performance  of  most  airborne  radar  systems  in  de¬ 
tecting  targets  on  or  near  the  wind-swept  surface  of  the  sea  is  return  echoes  from  the  sea 
surface  known  as  sea  clutter.  One  of  the  goals  in  this  study  is  to  develop  a  quantitative 
understanding  of  microwave  scattering  from  the  open  ocean  and  from  site-specific  coastal 
land  areas  for  medium-  and  high-PRF  airborne  surveillance  radar  systems  operating  in 
the  L-band  wavelengths.  Altitude  and  sidelobe  returns  from  sea  and  coastal  land  clutter 
play  important  roles  in  determining  the  minimum  signal-to-clutter  plus  noise  ratio  and 
clutter-to-noise  ratio  to  derived  the  required  antenna  sidelobe  levels  and  probability  of 
detection  [25].  It  is  therefore,  important  to  study  the  backscatter  from  the  sea  surface 
and  coastal  land  areas  for  a  wide-range  of  radar  grazing  angles,  in  particular  returns  from 
short-range,  near- normal  incidence  angle  [21].  This  report  evaluates  clutter  returns  from 
a  maritime  environment  for  airborne  surveillance  L-band  VV-polarised  radar  system  for  a 
wide  range  of  grazing  angles  (0.1°  <  4>  <  90°),  radar  look  directions,  varying  sea  conditions 
and  site-specific  coastal  land  areas. 

The  concepts  of  microwave  scattering  theory  of  slightly  rough  surface  have  been  applied 
to  study  the  backscattering  from  sea  surface  for  L-band  VV-polarisation  under  standard 
atmospheric  conditions  [16,  14].  This  allows  us  to  formulate  a  second-order  composite  sea 
clutter  model  or  MCSM  in  short,  to  predict  the  normalised  sea  backscattering  coefficient, 
a°,  for  varying  sea  conditions.  Based  on  the  two-scale  sea  surface  model  proposed  by 
Guinard  and  Daley  [16]  and  surface  scattering  model  proposed  by  Fung  and  Lee  [14],  the 
MCSM  sea  clutter  model  takes  into  account  polarisation,  wind  speed,  azimuth  and  grazing 
angles  variations.  In  this  model,  the  backscattering  coefficient,  cr°,  can  be  represented  as 
the  sum  of  separate  contributions  from  specular  return  and  resonant  Bragg-scattering. 

In  the  report,  we  also  briefly  review  three  commonly  used  semi-empirical  sea  clutter 
models.  They  are  the  Georgia  Institute  of  Technology  (GIT),  the  Hybrid  (HYB),  and 
the  Technology  Service  Corporation  (TSC)  sea  clutter  models.  These  sea  clutter  models 
are  collated  with  five  sets  of  measured  L-band  VV-polarised  data.  Of  the  five  data  sets 
used,  three  sets  of  calibrated  airborne  L-band  data  were  collected  by  the  Naval  Research 
Laboratory  (NRL)  in  the  late  60s  and  early  70s  from  several  sea  locations  and  in  varying 
sea  conditions.  The  fourth  set  of  data  was  compiled  by  Nathanson  and  the  updated  sea 
backscattered  coefficient  values  were  published  in  his  most  recent  book  [29].  The  last  set 
of  data  was  obtained  from  measured  L-band  VV-polarised  polarimetric  synthetic  aperture 
radar  (POLSAR)  images  of  the  Darwin  region.  These  images  were  acquired  during  the 
PACRIM  mission  by  the  Jet  Propulsion  Laboratory  on  the  23rd  November  1996. 

For  a  radar  system  operating  close  to  the  coastal  region,  its  detection  and  tracking 
capabilities  will  be  severely  affected  by  compounded  backscatter  from  the  coastal  water, 
foreshore  vegetation,  and  man-made  structures.  The  clutter  returns  from  the  coastal  areas 
will  depend  on  the  coastal  terrain  characteristics  e.g.,  vegetation  type,  presence  of  man¬ 
made  structures,  ground  surface  roughness  and  many  more  [18,  30].  The  availability  of 
the  Darwin  POLSAR  data  allows  us  to  study  the  characteristics  of  backscatter  signals 
from  coastal  water,  coastal  saline  flat  surfaces  and  mangrove  forest,  and  their  dependence 
on  grazing  angle  for  an  L-band  VV-polarised  radar  system.  Terrain  scatterer  types  in  the 
POLSAR  images  were  identified  using  the  TOPO-IOOK  and  TOPO-50K  maps. 
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2  Backscattering  from  the  Sea  Surface  in 
Standard  Atmospheric  Conditions  (cr°) 

The  magnitude  of  the  radar  signal  backscatter  from  the  sea  surface  is  dependent  on  a 
number  of  parameters;  these  include  the  wavelength  of  the  transmitted  microwave  signal, 
the  incidence  or  grazing  angle,  the  polarisation  of  the  transmitted  wave,  the  reflective  index 
of  the  sea  surface,  possible  surface  atmospherics  effect  such  as  surface  and  evaporation 
ducts,  and  environmental  conditions  (sea  state  and  wind  speed). 

2.1  Dependence  on  Radar  Grazing  Angle 

Studies  observing,  modelling  and  predicting  the  backscattering  coefficients,  cr°,  from 
sea  surface  are  probably  the  most  extensive  research  aspect  of  radar  sea  return.  One  of 
the  main  factors  affecting  signal  returns  from  the  sea  surface  is  the  grazing  angle  of  the 
transmitted  wave.  The  range  of  possible  grazing  angles  can  be  stratified  into  three  distinct 
regions  as  shown  in  Figure  2.1  [29,  25]: 


o 


Diffraction 

Trapping 


Figure  2.1:  Relationship  between  grazing  angle,  polarisation  and  radar  sea  backscatter 
returns 


1.  high- grazing- angle  region  (quasi-specular)  for  grazing  angle  between  ips  to  90°,  ip  6 
["0s  j  90°); 

In  the  high-grazing-angle  region,  the  magnitude  of  the  backscatter  signal  from  sea 
surface  is  controlled  by  mirror-like  reflection  of  the  local  flat  sea  surface  facets  tilted 
to  align  normally  to  the  incident  signal. 
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2.  plateau- grazing- angle  region  for  grazing  angle  between  critical  angle  to  ips,  ip  e 

tytM  ; 

The  magnitude  of  the  backscatter  signal  from  the  sea  surface  is  controlled  by  the 
Bragg  resonant  scatterers. 

3.  low-grazing-angle  region  for  grazing  angle  below  the  critical  angle,  ip  <  ipp, 

Breaking  waves,  wedge  scattering,  multipath  reflections  and  the  presence  of  ducting 
significantly  increase  the  backscatter  signal. 


The  critical  angle  is  [26] 


sin  ipt  = 


A 

4irhe 


(1) 


where  he  is  the  RMS  wave  height.  Note  that  the  boundaries  between  the  three  regions 
change  with  wavelength,  polarisation,  and  surface  condition  [22]. 


2.1.1  Quasi-Specular  Grazing  Angle  Region 

In  the  quasi-specular  region,  the  magnitude  of  the  backscatter  signal  is  controlled  by 
the  mirror-like  reflection  of  local  flat  facets  tilted  by  the  long  gravity  waves  [23,  50].  Local 
surface  waves  tilted  normal  to  the  radar  signal  cause  large  specular  returns  for  grazing 
angles  near  normal  incidence. 

<7°  =  (2) 

Specular  Returns 

The  main  contributor  to  the  specular  returns  are  the  long  gravity  waves  [13]. 


2.1.2  Plateau  or  Bragg  Backscattering  Region 


The  radar  signal  backscatter  from  the  sea  surface  for  grazing  angles,  ipt  <  ip  <  ips,  at 
L-band  is  predominantly  due  to  ocean  waves  with  wave  length  similar  to  the  transmitted 
wavelength.  Known  as  Bragg  scatterers  or  Bragg’s  resonant  ripples ,  these  water  waves 
(half  the  radar  wavelength)  are  the  main  source  of  backscattering  for  radar  operating  at 
grazing  angle,  ip  between  25°  and  75°  [42,  13,  44,  4],  The  relationship  between  Bragg’s 
resonant  sea  (or  spatial)  wavelength  and  transmitted  wavelength  is 


A  = 


A 

2  cos  ip 


(3) 


where  A  is  the  sea  wavelength,  A  is  the  transmitted  wavelength,  and  ip  is  the  local  grazing 
angle.  Lee  and  associates  [20]  studied  the  Bragg  scatterers  in  detail  and  concluded  that 
Bragg  scatter  is  a  single-bounce  return  in  the  absent  of  depolarisation  and  with  negligible 
cross  polarisation.  Note  that  the  Bragg  resonant  wavelength,  A,  is  shorter  at  smaller 
grazing  angle.  The  sea  wavenumber  is 


K  = 


2tt 

T 


(4) 
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The  first-order  Bragg  scattering  mechanism  for  L-band  radar  falls  in  the  short  gravity¬ 
capillary  waves  region  at  midrange  grazing  angles  in  the  absence  of  longer  and  higher 
gravity  waves. 

These  short  gravity-capillary  waves  provide  the  mechanism  for  microwave  backscatter- 
ing  by  forming  distributed  facets  atop  the  gently  undulating  long  gravity  waves  to  reflect 
the  signal  back  to  the  radar  system.  The  long  gravity  waves  have  four  significant  effects 
on  the  short  waves:  (1)  tilting  the  scattering  surface,  thereby,  modifying  the  local  grazing 
angle  [16,  38]  as  shown  in  Figure  2.2;  (2)  local  straining,  causing  amplitude  variation  over 


Short  Gravity  Waves 


Figure  2.2:  Two-scale  Sea  Surface 

the  longer  wave,  with  maximum  straining  occurring  near  the  crest  of  the  gravity  wave  [38]; 
(3)  advection  by  the  orbital  velocity  of  the  gravity  wave  [38];  and  (4)  variation  in  radar 
range  phase  path  [38]. 


2.1.3  Low  Grazing  Angle  Region 

At  small  grazing  angles,  the  sea  surface  appears  in  the  electromagnetic  sense  to  be 
smooth  for  an  L-band  radar.  There  are  a  number  of  scattering  mechanisms  that  contribute 
to  the  radar  return  signal  from  the  sea  surface  at  low  grazing  angles,  in  particular,  breaking 
waves  due  to  the  large  gravity  waves,  the  effects  of  ducting,  multipath  reflection,  and  wedge 
scattering. 

“Wedge  scattering”  is  the  scattering  contribution  from  sea  surface  with  small  radius 
of  curvature  relative  to  the  radar  wavelength  [24].  Lyzenga  and  associates  [24]  proposed 
an  analytic  expression  to  consider  the  contribution  of  wedge  scattering  to  the  sea  surface 
backscattering  coefficient,  a°. 


o  _  92ln 

t^VV -wedge  ~~  Z 


(5) 


where  N  is  the  number  of  wedge  crests  per  unit  area  and  L  is  their  average  length. 
However,  it  has  also  been  noted  by  the  authors  that  they  are  unsure  how  these  quantities 
can  be  measured  accurately. 
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2.2  Dependence  on  Wind  Speed 

The  <7°  variation  of  the  sea  surface  is  dependent  on  the  wind  speed.  The  relation  is 
however  complex  and  uncertain.  In  early  publications,  Valenzuela  and  associates  [46]  have 
noted  that  the  wind  dependence  of  the  sea  spectrum  is  identical  to  the  wind  dependence  of 
the  a°  of  the  sea.  Using  the  NRL-4FR  data,  Valenzuela  [46]  concluded  that  the  a°  is  less 
sensitive  at  L-band  to  increasing  wind  speed  and  the  relationship  is  best  approximated  by 
a  power  law.  Similarly,  Daley  [7]  and  Wetzel  [39]  have  noted  that  the  relationship  between 
<7°  from  sea  and  wind  speed  for  L-band  is  given  by  a  power-law  of  the  form 

a°  «  U2v  (6) 

where  the  estimated  values  of  2 v  from  the  NRL-4FR  data  [7]  is 


0.4, 

Upwind 

0.2, 

Downwind 

(7) 

-0.1, 

Crosswind 

Similarly,  Daley  [9]  has  reported  that  the  measured  a°  (JOSS  I  data)  is  consistent  with  the 
power  law  described  by  Valenzuela  [46].  Thompson  and  associates  [41]  have  also  reported 
the  wind-speed  exponent  to  be  approximately,  2u  =  0.5.  These  data  indicate  a  dependence 
upon  surface  wind  speed  since  the  L-band  backscatter  signal  is  approximately  proportional 
to  the  square  root  of  the  wind  speed. 

However  Donelan  and  Pierson  [13]  have  argued  that  the  power-law  relationship  un¬ 
derestimates  the  backscatter  for  high  wind  speed  and  overestimates  the  backscatter  at 
low  wind  speed.  They  have  shown  that  the  high-wave  number  gravity-capillary  waves  is 
a  function  of  both  the  wind  speed  and  viscosity  (variation  of  viscosity  of  sea  water  with 
temperature  and  salinity). 

In  addition,  they  have  also  proposed  a  threshold  wind  speed,  whereby  no  backscattering 
in  the  plateau  region  (mid-range  incidence  angles)  is  observed  below  the  threshold,  and 
that  this  threshold  speed  is  a  function  of  the  incidence  angle  and  the  sea  water  temperature. 
This  minimum  wind  speed  is  required  to  overcome  viscosity  damping  and  allow  the  short 
gravity-capillary  waves  to  grow  [34].  At  L-band  for  30°  C  sea  water  temperature  and 
ip  —  70°,  Bragg  resonant  backscattering  is  observed  when  the  wind  exceeds  1.7  m/s.  At 
ip  =  25°,  it  only  need  to  exceeds  1.61  m/s  [13]. 

The  data  reported  by  Morris  [26]  in  Figure  2.3  shows  that  for  VV  polarised  sea  clut¬ 
ter,  there  appears  to  be  little  dependence  on  wind  speed  for  grazing  angles  in  the  low 
and  plateau  regions.  However,  a  slight  dependence  on  wind  speed  is  observed  from  HH 
polarised  sea  clutter  at  low  grazing  angles. 

We  can  conclude  that  in  L-band  (VV  polarisation),  the  backscatter  from  the  sea  surface 
has  the  following  behaviour: 

1.  The  radar  backscatter,  <7°,  from  sea  increases  slightly  with  increasing  wind  speed. 

2.  Wind  dependence  increases  with  decreasing  grazing  angle,  and  at  near  normal  inci¬ 
dence,  the  backscatter  is  inversely  proportional  to  wind  speed  [14,  19]. 
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Relationship  between  Backscattered  Coefficient,  Polarisation,  Sea  State  and  Wind  Speed 


2.3  Dependence  on  Wind  Direction 

The  results  published  from  several  trials  have  shown  that  the  backscattered  signal  from 
the  sea  surface  is  dependent  on  a  factor  that  combines  both  the  wind  magnitude  with  its 
relative  direction  to  the  boresight.  The  sea  clutter  is  strongest  when  viewed  in  the  upwind 
direction,  intermediate  in  the  downwind  direction,  and  weakest  in  the  crosswind  direction 
for  S,  X,  and  Ku-band  sea  clutter  [39]. 

However,  for  an  L-band  radar,  the  sea  surface  appears  smooth,  resulting  in  negligible 
difference  between  upwind,  and  downwind  backscatter  returns  [14].  Similar  findings  were 
reported  by  Chan  [4]  and  Thompson  and  associates  [41]  for  L-band  radar.  Chan  found  that 
the  upwind/downwind  ratio  of  a°  was  close  to  unity,  and  the  is  very  high  for  radar 
operating  within  these  frequency  bands.  Thompson  and  associates  [41]  found  little  or  no 
dependence  of  L-band  backscatter  on  wind  direction  for  stable  atmospheric  conditions,  but 
observed  that  a  strong  dependence  on  wind  direction  occurred  during  unstable  conditions. 


2.4  Dependence  on  Rain 

Rain  does  not  produce  an  overriding  influence  with  measurement  of  backscatter  char¬ 
acteristics  of  sea  surface  for  L-band  radar  [50].  In  addition,  studies  have  shown  that  the 
VV-polarised  sea  surface  backscatter  is  insensitive  to  rain  rate  [39] . 
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3  Data  Source  for  L  Band 

Five  sets  of  measured  L-band  VV-polarised  data  have  been  used  to  test  the  accu¬ 
racy  of  the  GIT  ,  HYB  ,  TSC  and  the  proposed  MCSM  sea  clutter  models  in  predicting 
sea  backscattering  values  for  varying  sea  conditions  and  grazing  angles.  Measured  sea 
backscattered  signals  on  which  analyses  in  this  report  are  based  on  are:  (1)  Naval  Re¬ 
search  Laboratory  (NRL)  Four-Frequency  radar  (4FR)  [46,  9]  calibrated  measurements; 
(2)  measurements  compiled  by  Nathanson  [29];  and  (3)  Jet  Propulsion  Laboratory  POL- 
SAR  Darwin  data.  Brief  descriptions  of  the  NRL  4FR  sets  of  data,  Nathanson’s  data  and 
the  DARWIN  POLSAR  data  are  summarised  below. 


3.1  Naval  Research  Laboratory  Airborne  Trials  [46,  16,  8,  9] 

Measurements  of  the  normalised  median  radar  backscattering  coefficient,  cr^,  from 
the  sea  surface  for  varying  grazing  angles  and  sea  conditions  were  recorded  by  the  Naval 
Research  Laboratory  (NRL)  in  the  60s  and  early  70s  [39,  8,  9].  The  NRL  deployment 
utilised  an  airborne  system  equipped  with  a  four-frequency  (4FR)  radar  system,  trans¬ 
mitting  a  sequence  of  frequencies  alternately  in  horizontal  and  vertical  polarisation.  The 
carrier  frequencies  were:  (1)  UHF  (428  MHz);  (2)  L-band  (1.228  GHz);  (3)  C-band  (4.455 
GHz),  and  (4)  X-band  (8.910  GHz).  The  system  was  calibrated  against  standard  metal 
spheres,  and  measurements  were  recorded  for  upwind,  downwind,  and  crosswind  in  wind 
speeds  ranging  from  5  to  50  knots  (up  to  7  -8  m  wave  height),  and  over  grazing  angles 
ranging  from  5°  to  90°.  Corresponding  wind  speeds  and  wave  heights  were  recorded  in 
the  measurements  areas  from  instrumented  ships  [39,  8,  9]. 

1.  NRL-4FR  Data  [46] 

Sea  backscattered  measurements  recorded  using  the  NRL  4FR  radar  system,  off  New 
Jersey  in  December  1964,  off  Puerto  Rico  in  July  1965  and  in  the  North  Atlantic 
in  February  1969  were  considered  in  this  report.  Surface  truth  information  for  the 
Puerto  Rico  deployment  was  measured  by  a  team  of  researchers  from  the  Applied 
Physics  Laboratory,  John  Hopkins  University,  in  an  instrumented  vessel  [16].  For  the 
North  Atlantic  deployment,  the  surface  truth  information  was  provided  by  oceano¬ 
graphic  surface  vessels  operated  by  the  French  and  British  governments  [16].  We 
will  hereafter  refer  the  data  measured  as  the  NRL-4FR  data. 

2.  Joint  Ocean  Surface  Study,  1970  (JOSS  I)  [9] 

The  NRL  4FR  system  was  employed  by  the  Naval  Oceanographic  Office  in  Jan¬ 
uary  1970  to  measure  calibrated  sea  backscattered  signals  for  various  sea  states, 
radar  wavelength,  polarisation,  and  depression  angle  off  the  vicinity  of  Argus  Island, 
Bermuda. 

The  4FR  system  was  operating  at  a  pulse  width  of  0.5/u.s,  PRF  of  683  Hz,  IF  band¬ 
width  of  10  MHz,  and  range-gate  width  of  24  ns.  Surface  truth  information  was 
acquired  utilising  an  altimeter  and  aerial  cameras  [9].  Wind  and  wave  height  infor¬ 
mation  were  recorded  by  observers  stationed  on  Argus  Island. 
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3.  Joint  Ocean  Surface  Study,  1971  (JOSS  II)  [10] 

The  NRL  4FR  system  was  employed  by  the  Naval  Oceanographic  Office  in  February 
1971  to  measure  calibrated  sea  backscattered  signals  for  various  sea  states,  radar 
wavelength,  polarisation,  and  depression  angle  off  the  East  Coast  of  Northern  Amer¬ 
ica. 

Surface  truth  data  was  acquired  at  the  site  of  study  from  an  instrumented  buoy 
and  an  ocean  station  vessel.  Wind  velocity,  direction,  and  average  wave  height 
information  were  obtained  hourly  from  both  of  these  instruments  [10].  The  values  of 
the  4FR  system  used  during  the  JOSS  II  program  were  PRF,  512  Hz;  pulse  width, 
1.0  /zs;  IF  bandwidth,  10  MHz;  and  range-gate  width,  24  ns  [10]. 

For  all  analyses  carried  out  in  this  report,  the  median  sea  backscattered  values  were 
converted  to  mean  values  assuming  an  exponentially  distribution  or  Rayleigh  in  amplitude 

[39]- 

3.2  Normalised  Mean  Backscattering  Coefficient  (a0)  com¬ 
piled  by  Nathanson  [29] 

Nathanson  compiled  measurements  of  the  sea  surface  normalised  mean  backscattering 
coefficient,  cr°,  given  the  sea  state,  polarisation,  and  carrier  frequency  from  approximately 
60  experiments.  Some  of  the  values  tabulated  in  his  most  recent  book  [29]  were  approx¬ 
imated  by  extrapolation  and  interpolation  between  points  or  were  averaged  of  upwind, 
downwind,  and  crosswind  situations.  Grazing  angles  were  adjusted  from  depression  or 
incidence  angles,  and  in  some  data,  median  values  were  adjusted  to  mean  values. 

Some  general  comments  about  the  database: 

(a)  Measurements  were  collected  from  L-band  radar  systems  operating  at  pulse  widths 
ranging  from  0.5-  to  5-^zs. 

(b)  The  database  contains  many  uncertainties  at  low  carrier  frequencies  and  low  grazing 
angles. 

(c)  Nathanson  believed  that  the  a°  values  reported  for  low  grazing  angles  to  be  higher 
than  true  mean  values  due  to  the  following: 

(i)  Early  data  were  generated  from  less  sensitive  radar  and  values  reported  were 
only  of  data  above  normal  (above  known  mean  values),  and  from  unknown  wind 
speed,  sea  state,  and  ducting  conditions; 

(ii)  The  high  occurrence  of  ducting  increases  the  a°  values  at  very  low  grazing 
angles. 

(d)  Although  the  mean  values  published  in  his  most  recent  book  represent  lower  a° 
values  than  that  reported  in  1969,  they  might  still  be  higher  than  true  mean  values 
due  to  ducting  conditions. 

Although  the  data  compiled  by  Nathanson  is  biased  above  the  true  mean  values,  it 
nevertheless  provides  the  upper  limit  scenario  when  evaluating  the  radar  performance  in 
average  sea  conditions. 
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3.3  POLSAR  Darwin  Data 

POLSAR  is  a  NASA-sponsored  polarimetric  synthetic  aperture  radar  system  developed 
at  the  Jet  Propulsion  Laboratory  in  California.  The  POLSAR  system  provides  high  quality 
polarimetric  SAR  data  in  three  frequencies;  P  (A  =  67.0  cm,  /  =  0.44  GHz),  L  (A  =  23.9 
cm,  /  =  1.25  GHz)  and  C  (A  =  5.6  cm,  /  =  5.39  GHz)  bands  in  four  transmit /receive 
polarisation  (VV,  HH,  VH,  HV). 

3.3.1  Study  Area 


Figure  3.1:  Location  of  the  Study  Area 

The  study  area  of  interest  was  conducted  in  the  Darwin  region  (130°  45  ’  B,  12  °  30’  S) 
located  in  the  Northern  Territory  of  Australia  (Figure  3.1).  The  monsoon  season  in  this 
region  falls  between  December  and  March,  with  average  rainfall  of  approximately  1624 
mm  [3].  The  daily  maximum  temperature  ranges  between  30  °C  to  34°C  (hottest  months 
are  November  and  December),  and  the  minimum  temperature  ranges  between  19° C  to 
27°  C  (coolest  month  is  July)  [3]. 

The  study  area  consists  of  several  land  cover  classes,  but  only  vegetation  close  to  the 
coastal  region  is  considered  in  this  report.  The  Northern  Territory  coast  can  be  seen  as  an 
heterogeneous  evergreen  mangrove  forest  covering  an  area  of  4,120  km2  [3].  The  mangrove 
forest  of  the  Northern  Territory  forms  a  distinct  lineal  growth  along  the  shorelines  of  the 
coastal  regions.  Mangrove  community  patterns  surrounding  the  coastal  areas  of  Darwin, 
with  a  total  extent  of  20,4000  hectares  is  dominated  by  Ceriops  tagal  forest  as  shown 
in  Figure  3.3  (a).  The  height  of  these  trees  ranges  from  1  to  10  meters  and  very  often 
homogeneous  population  of  this  tree  type  can  be  found  in  the  mid  and  high  tidal  flat 
regions  and  in  the  hinterland  regions,  covering  40%  of  the  mangrove  forest  areas  (Figure 
3.2)  .  Soil  composition  in  these  regions  are  generally  muds  and  sandy  muds. 
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(a)  (b) 

Figure  3.3:  (a)  Ceriops  tagal  low  closed-forest  ( mid  tidal  flat),  and  (b)  Rhizophora  stylosa 
closed- forest/low  closed-forest  (shoreline  forest)  [3] 


Rhizophora  stylosa  closed-forest  is  also  commonly  found  along  the  tidal  creeks  and  the 
seaward  fringes  of  the  main  tidal  flat  [3],  varying  in  height  from  6  to  18  meters  (covering 
approximately  33%  of  the  mangrove  areas).  As  shown  in  Figure  3.3  (b),  these  trees 
generally  have  straight  single  bole  in  the  tallest  stumps  with  aerial  root  structure  forming 
dense  entanglement  visible  during  low-tide.  The  height  of  these  mangrove  trees  decreases 
as  we  move  towards  a  higher  tidal  level.  It  has  been  noted  that  the  mangrove  vegetation 
of  the  above  two  species  are  particularly  sensitive  to  cyclonic  activities  common  in  the 
Northern  Territory.  Defoliation  and  high  mortality  is  observed  following  this  disturbance 
[3] ,  resulting  in  a  different  signature  of  backscatter  characteristics  observed  by  the  radar 
system.  Mud  flats  on  the  seaward  margin  are  generally  devoid  of  vegetation  below  the  3 
m  tidal  level  [3]. 

The  JPL  POLSAR  system  was  flown  over  the  study  site  onboard  a  DC-8  on  the  23rd 
November  1996.  Two  flight  passes  were  selected  and  the  two  L-band  polarimetric  SAR 
images  of  vertical  (VV)  polarisation  collected  with  bandwidth  of  40  MHz  are  shown  in 
Figure  3.4  and  3.5. 

The  POLSAR  data  were  processed  as  a  60km  strip  image  by  the  JPL’s  Radar  Data 
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Center  using  version  6.10  (CM5512  -  Figure  3.4)  and  5.10  (CM5148  -Figure  3.5)  integrated 
processor.  The  number  of  looks  processed  in  the  azimuth  is  18  and  1  in  range. 


Table  3.1:  POLSAR  Specifications  for  the  23  November  1996  Darwin  Deployment 


Description 

CM5512 

CM5148 

Acquisition  Time  (GMT) 

17378.0 

16746.0 

Frequency  (GHz) 

1.218 

1.218 

Near  Slant  Range  (m) 

9013.26 

9013.26 

Far  Slant  Range  (m) 

17537.36 

19244.51 

Near  Look  Angle  (9i) 

22.5° 

22.6° 

Far  Look  Angle  ( 9i ) 

61.7° 

61.7° 

Range  Resolution  (m) 

3.3310 

3.3310 

Azimuth  Resolution  (m) 

9.2592 

9.2592 

Image  Size  (pixels) 

4183  x  2560 

1095  x  2560 

Limited  ground  truth  is  available  for  vegetation  measurements  of  the  site  of  interest. 
However,  detail  ground  truthing  of  moored  vessels,  sea  and  meteorological  conditions  were 
conducted  by  a  team  of  DSTO  scientists.  On  the  23rd  November  1996,  moored  vessels 
were  photographed  and  position  recorded  at  the  Fannie  Bay  (by  Mr  David  Clark)  and 
Port  Darwin  (by  Dr.  Ed  Kruzins)  areas.  At  the  same  time,  aerial  photographs  of  the 
POLSAR  swath  area  were  obtained  from  an  altitude  of  26,000  ft.  with  the  assistance  of 
the  local  coastwatch  personnel.  Two  passes  of  the  Darwin  harbour  were  conducted  by  Mr 
Kim  Meaney:  (1)  an  approximately  east-west  pass;  and  (2)  an  approximately  north-south 
pass.  Long  focal  lenses  of  600mm  and  1200mm  were  used  to  ground  truth  moored  vessels 
in  the  Darwin  harbour  (look  angle  to  targets  is  approximately  45°).  In  addition,  a  4  ft. 
square  corner  reflector  was  placed  in  an  open  area  at  the  northern  end  of  the  Darwin 
airport.  Recorded  measurements  of  the  sea  and  meteorological  conditions  are  tabulated 
in  Table  3.2. 


Table  3.2:  Hydro-meteor  Parameters  for  the  23  November  1996  Darwin  Deployment  (Pro¬ 
vided  by  Dr.  Ed  Kruzins,  DSTO) 


Description 

Darwin  Harbour 

Fannie  Bay 

Acquisition  Time  (CST) 

12:17  pm 

10:40  am 

Air  Temperature  (Dry) 

35.0°C 

32.5°C 

Air  Temperature  (Wet) 

29.0°C 

27.0°C 

Relative  Humidity 

68% 

64% 

Cloud  Cover 

5/8  Octa 

3/8  Octa 

Cloud  Type 

Cirrus/Cumulus 

Cirrus  and  Strato-Cumulus 

Wind  Speed 

14  km/h 

4  kt 

Wind  Direction 

170° 

SE 

Wave  Height 

1  cm 

Calm 
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Figure  3.4:  POLSAR  Darwin  L-band  VV  data  (CM5512) 
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Figure  3.5:  POLSAR  Darwin  L-band  VV  data  (CM5148) 
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3.3.2  SAR  Terrain  Analysis  Program 

A  program  has  been  written  in  MATLAB  to  analyse  the  JPL  POLSAR  images.  The 
main  user  interface  is  shown  in  Figure  3.6. 


Figure  3.6:  SAR  Analysis  GUI  Interface 

A  number  of  flexible  options  have  been  incorporated  into  the  software  design.  These 
options  allow  the  user  to  specify  the  near  and  far  slant  ranges,  radar  altitude,  and  near 
and  far  field  directions.  In  addition,  for  high  resolution  data  or  images,  the  user  is  able  to 
view  the  image  in  sub-quadrant  (maximum  of  4) . 

The  program  allows  the  user  to  operate  in  two  modes: 

1.  Pixel  Mode 
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In  the  pixel  mode,  a  dialog  box  pops-up  allowing  the  user  to  view  the  backscattering 
coefficients,  <r°,  of  individual  pixel  and  grazing  angle  at  the  mouse  position.  An 
example  of  the  dialog  box  is  shown  below  (Figure  3.7). 


Figure  3.7:  SAR  Pixel  Analysis  GUI  Interface 

2.  Area  Mode 

(a)  Segment  Profile 

The  segment  profile  mode  allows  the  user  to  specify  a  line  segment  of  the 
image  of  interest  using  a  mouse.  The  program  calculates  the  grazing  angle  and 
normalised  backscattering  coefficient,  a°,  for  every  pixel  on  the  line  segment. 

(b)  Area  Profile 

The  area  profile  mode  allows  the  user  to  crop  a  small  area  of  the  image  to 
determine  a  set  of  backscattering  coefficient,  <r°,  with  its  corresponding  grazing 
angle  value.  Provisions  have  also  been  made  to  allow  the  user  to  view  the 
histogram  of  a°  of  the  selected  area,  and  save  the  measured  values  to  a  file. 

Radar  backscattered  signals  from  the  sea  surface  were  extracted  from  the  CM5512 
and  CM5148  POLSAR  images  around  the  Port  Darwin  area  at  a  distance  away  from  the 
coastal  saline  flat  surfaces.  Figure  3.8  shows  the  o°  as  a  function  of  grazing  angle. 


I 

I 


Figure  3.S:  Normalised  a°  coefficients  from  the  sea  surface  (L-VV,  U  =  3.88  m/s) 


Darwin  Su  Raflacthrlty  Modal  (L-band.  W) 


Grazing  Anglo  (degrooa) 
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4  Semi-Empirical  Sea  Clutter  Models 

There  are  a  number  of  established  semi-empirical  sea  clutter  models  available  for  de¬ 
termining  the  sea  backscattering  coefficient,  a°,  given  knowledge  of  the  radar  parameters, 
sea  conditions,  and  method  of  observation.  These  models  provide  a  fairly  accurate  estima¬ 
tion  of  the  sea  clutter  backscattering  coefficient,  o° .  Three  empirical  sea  clutter  models 
suitable  for  determining  the  radar  sea  returns  for  L-band  radar  have  been  selected  for 
further  consideration  in  this  report.  The  first  model  is  derived  from  the  Georgia  Institute 
of  Technology  (GIT),  the  second  model  is  a  Hybrid  (HYB)  model  that  include  works  from 
Barton  [35],  and  the  third  model  is  derived  from  Technology  Service  Corporation  (TSC) 
[1],  These  three  sea  clutter  models  are  applicable  to  a  wide  range  of  frequencies.  However, 
the  GIT  (0.1°  <  ip  <  10°)  and  HYB  (0.1°  <  ip  <  30°)  models  are  applicable  only  to 
low-grazing  angles,  while  the  TSC  (0.1°  <  <  90°)  model  is  applicable  to  a  wide  range 

of  grazing  angles  [1]. 

The  effectiveness  of  the  GIT  ,  HYB  ,  and  the  TSC  models  in  modelling  sea  clutter 
were  compared  with  measured  data  under  similar  sea  conditions.  Radar  sea  clutter  mea¬ 
surements  are  derived  from:  (1)  NRL-4FR  data  [46,  8],  and  (2)  a°  measurements  compiled 
by  Nathanson  [29]. 

The  GIT  ,  HYB  and  TSC  models  are  described  below: 

(a)  Georgia  Institute  of  Technology  (GIT)  Sea  Clutter  Model  [37,  1] 

The  Georgia  Institute  of  Technology  Sea  Clutter  Model,  hereinafter  referred  to  as 
the  GIT  clutter  model,  is  derived  from  multiple  data  sets  observed  under  “stan¬ 
dard”  propagation  conditions, and  mathematical  models  of  the  various  sea  surface 
backscattering  mechanism.  The  GIT  sea  clutter  model  describes  the  backscattering 
coefficient  as  a  function  of  three  factors:  (1)  Interference  factor,  Gm:  (2)  sea  condi¬ 
tion  -  wind  speed  factor,  Gu;  and  (3)  direction  dependence  -  wind  direction  factor, 
Gw.  The  interference  factor  uses  wave  height  to  take  into  account  the  transition 
in  receiver  power  from  R~ 3  (plateau  region)  dependency  to  R~7  in  the  low-grazing 
angles  [37].  This  model  is  listed  in  Appendix  B.l. 

(b)  Hybrid  (HYB)  Sea  Clutter  Model  [35] 

The  hybrid  model  (HYB)  is  derived  from  information  of  the  average  sea  clutter  data 
compiled  by  Nathanson  in  1969,  works  by  Barton,  and  features  of  the  GIT  model 
[29].  The  hybrid  model  is  suitable  for  assessing  the  sea  backscattered  coefficient,  a° , 
for  grazing  angle  between  0.1°  and  30°.  This  model,  hereinafter  referred  to  as  the 
HYB  model,  is  listed  in  Appendix  B.2. 

(c)  Technology  Service  Corporation  (TSC)  Sea  Clutter  Model  [1] 

The  Technology  Service  Corporation  sea  clutter  model,  hereinafter  referred  to  as  the 
TSC  sea  clutter  model,  is  an  empirical  model  derived  from  Nathanson’s  average  sea 
clutter  data  published  in  1969.  The  TSC  model  assumes  that  the  Nathanson’s  data 
is  the  average  for  all  look  directions,  and  hence  treats  it  as  cross-wind  data.  The 
TSC  models  takes  into  consideration  the  sea  conditions  such  as  Douglas  sea  states 
to  compute  wind  speed  and  wave  height. 
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It  has  been  noted  by  Reilly  [35]  that  the  TSC  model  includes  the  effect  of  some 
ducting  and  hence  predicts  higher  a°  at  low-grazing  angles.  Dockery  [12]  suggested 
that  the  in  absence  of  prior  information  about  propagation  conditions,  the  TSC 
model  better  represents  average  sea  conditions  [37,  12].  This  model  is  listed  in 
Appendix  B.3. 

A  summary  of  the  semi-empirical  sea  clutter  capabilities  is  tabulated  below  [37,  1]: 


Table  4-1:  Semi- Empirical  Sea  Clutter  Capabilities 


Description 

GIT 

HYB 

TSC 

Frequency  Band  (GHz) 

1-100 

0.5  -  35 

0.5  -  35 

Sea  Conditions 

Wind 

Douglas  Sea  States 

Douglas  Sea  States 

Geometry 

Speed  (m/s) 

(0-5) 

(0-5) 

Grazing  Angle  (°) 

0.1  -  10 

0.1  -  30 

0.1  -  90 

Aspect  Angle  (°) 

0  -  180 

0  -  180 

0-  180 

Polarisation 

VV,  HH 

VV,  HH 

VV,  HH 

4.1  The  Use  of  Nathanson  and  NRL-4FR  Data  Sets  to  Test 
the  GIT,  HYB,  and  TSC  Models 

The  average  backscattering  values  compiled  by  Nathanson  and  the  NRL-4FR  data  sets 
have  been  used  to  test  the  accuracy  of  the  GIT  ,  HYB  ,  and  TSC  sea  clutter  models.  The 
plots  of  Figure  4.1  (a)  and  (b)  compare  these  three  sea  clutter  models  at  wind  speed  of  5.6 
m/s  and  12.1  m/s.  These  plots  apply  to  an  L-band  radar  with  carrier  frequency  of  1.25 
GHz  with  radar  transmitting  and  receiving  at  vertical  polarisation  and  an  upwind  radar 
azimuth  direction. 

(a)  Discussion  on  GIT  Sea  Clutter  Model 

The  GIT  model  takes  into  account  the  multipath  returns  at  low  grazing  angle.  The 
multipath  or  interference  factor,  Gm,  is  a  theoretically  derived  factor  to  account 
for  multipath  interference  for  Gaussian  distributed  wave  heights  [32].  Therefore,  at 
small  grazing  angle,  the  GIT  model,  exhibits  a  rapid  falloff  at  the  rate  of  ip4  [35] 
as  shown  in  Figure  4.1  (a)  and  (b).  This  small  angle  region  is  also  known  as  the 
“Interference  Region”.  However,  this  approach  is  justified  only  in  the  absence  of 
ducting  or  other  non-standard  propagation  effects  [12].  The  remaining  two  factors, 
Gu  and  Gw,  were  derived  empirically  to  account  for  wind  speed  and  wind  direction 
variation,  respectively. 

When  comparing  the  predicted  a°  with  the  measured  NRL-4FR  and  Nathanson’s 
data,  the  GIT  model  has  close  agreement  with  both  sets  of  data  for  ip  >  1°.  Since  the 
NRL-4FR  data  is  limited  to  grazing  angle  larger  than  10°,  it  is  difficult  to  evaluate 
the  accuracy  of  the  GIT  model.  However  for  grazing  angle,  ip  <  1°,  close  agreement 
between  Nathanson  data  and  GIT  predictions  occurs  only  at  higher  wind  speed 
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(12.1  m/s)  [refer  to  Figure  4.1  (b)].  Poor  agreement  can  be  seen  in  Figure  4.1  (a)  for 
low  surface  wind  speed.  Large  differences  between  the  GIT  model  and  Nathanson 
data  may  be  attributed  to  the  heavy  averaging  done  by  Nathanson  [49,  32]  on  data 
measured  under  a  wide  variety  of  propagation  conditions,  and  the  selective  recording 
of  measured  data. 

One  major  drawback  of  the  GIT  models,  is  the  limited  range  of  grazing  angles 
(0.1°  <  ip  <  10°).  Since  we  are  interested  in  determining  the  backscattered  returns 
from  sea  clutter  for  grazing  angles  between  ip  =  0°  and  ip  =  90°,  the  GIT  model  is 
not  suitable  for  our  application. 

(b)  Discussion  on  TSC  and  HYB  Sea  Clutter  Models 

The  TSC  and  HYB  sea  clutter  models  were  derived  by  empirically  fitting  curves  to 
large  quantities  of  Nathanson’s  data,  measured  under  a  wide  variety  of  propagation 
conditions  [12].  The  empirical  parameters  in  the  TSC  and  HYB  models  were  ad¬ 
justed  using  measured  a°  to  reduce  the  falloff  at  small  grazing  angles  [35].  These 
models  exhibit  a  ip2  falloff  behaviour  at  low  grazing  angles  [35].  As  shown  in  Fig¬ 
ure  4.1  (a)  and  (b),  the  sea  backscattering  coefficient,  <x°,  predicted  by  the  TSC 
and  HYB  model  is  significantly  higher  than  the  corresponding  GIT  predictions  for 
grazing  angles  less  than  1°.  These  sea  clutter  models  are  in  reasonable  agreement, 
especially  with  Nathanson  data  for  low  wind  speed  (5.6  m/s),  but  deviate  for  higher 
wind  speed  (12.1  m/s)  for  ip  <  1°. 

At  higher  grazing  angles  (ip  >  1°),  differences  between  these  models  and  Nathansons’ 
data  are  significantly  reduced.  In  the  absence  of  sea  ducting  conditions,  Dockery 
[12]  suggested  using  the  HYB  or  TSC  model  to  represent  the  average  propagation 
conditions.  However,  if  propagation  information  is  available,  Dockery  [12]  suggested 
using  a  modified  GIT  model  that  explicitly  includes  the  propagation  factor  (ducting 
or  nonstandard  atmospheric  effects)  to  model  the  sea  clutter.  Using  sea  clutter 
measurements  collected  with  a  2.9  GHz  radar,  upwind,  and  wind  speed  of  6.2  m/s, 
Dockery  [12]  showed  that  the  modified  GIT  model  is  in  reasonable  agreement  with 
the  measured  cr°  compared  with  the  HYB  or  the  GIT  models. 

Poor  agreement  between  the  measured  NRL-4FR  data,  and  the  predicted  a°  values 
were  observed  in  Figures  4.1  (a)  and  (b)  for  higher  grazing  angles  (ip  >  10°)  for  both 
wind  speeds.  Both  sea  clutter  models  underestimated  the  predicted  backscattered 
coefficient,  a°,  for  both  wind  speeds.  The  discrepancy  between  the  measured  and 
predicted  sea  backscattered  coefficients  increases  with  increasing  grazing  angles  and 
wind  speeds.  However,  in  the  quasi-specular  region  (ip  >  60°),  there  is  good  agree¬ 
ment  between  measured  backscatter  values  and  predicted  TSC  sea  clutter  values. 
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5  Composite  Sea  Clutter  Models 

The  composite  sea  clutter  models  to  be  developed  in  the  later  sections  are  based  upon 
the  two-scale  sea  clutter  model  developed  by  Guinard  and  Dailey  [16].  These  models  take 
into  account  polarisation,  wind  speed,  and  azimuth  plus  grazing  angles  variation. 

The  backscattering  coefficient,  a°,  from  the  sea  surface  can  be  represented  as  the  sum 
of  separate  contributions  from  specular  return  and  Bragg-scattering. 

a°  =  05  +  o' g  (8) 

Specular  Returns  Bragg  Scattering 

The  composite  sea  clutter  model  will  be  constructed  based  on  these  two  forms  of  scattering 
mechanism.  The  primary  scattering  mechanism  in  the  low-grazing  angle  and  plateau 
regions  is  assumed  to  be  Bragg  scattering  whereby  the  backscattering  coefficient,  is  in 
proportion  to  the  spectral  density  of  the  short  gravity-capillary  waves  for  a  fully  developed 
sea  [13].  At  near  normal  incidence,  the  effects  of  specular  return,  o°s,  due  to  the  stationary 
facet  normal  to  the  radar  signal  is  included. 


5.1  Specular  Backscattering  for  Fully  Developed  Sea 

Near-normal  incidence  backscatter  from  the  sea  surface  is  better  described  by  the 
specular  point  or  tilted  facet  model  [44,  47].  The  radar  backscatter  returns  from  the  sea 
surface  at  near  normal  incidence  is  in  proportion  to  the  number  of  mini-facets  normal  to 
the  incident  ray.  According  to  Barrick  [2],  specular  return  from  the  sea  surface  is  valid 
only  for  grazing  angles  greater  than  75°.  For  grazing  angles  less  than  75°,  the  backscatter 
return  is  more  diffuse  and  is  best  described  by  the  Bragg  diffraction  model. 

For  an  isotropic  rough  surface  of  Gaussian  statistics,  the  normalised  backscatter  return 
at  near  normal  incidence  is  [44] 

a%  =  sec2  exp(-  tan2  di/s 2)  (9) 

s 1 

where  i?(0)  is  the  Fresnel  reflection  coefficient  of  the  sea  surface  at  normal  incidence,  .s2 
is  the  mean-square  sea  slope  and  Oi  is  the  angle  of  incidence. 

Tilting  of  the  sea  surface  due  to  the  presence  of  swells  or  large  gravity  waves  would 
cause  the  actual  or  local  grazing  angles  to  be  different  from  those  calculated  from  a  flat 
sea  surface  assumption.  Under  the  flat  sea  surface  assumption,  the  relationship  between 
the  incidence  angle  and  the  grazing  angle  is 

=  90°  -  di  (10) 

If  the  sea  surface  is  elevated,  using  the  equivalent  earth  radius  method,  the  local  grazing 
angle  as  shown  in  Figure  5.1  is 

ip'  —  ip  +  dr  (11) 

where  dr  is  slope  of  the  sea  surface  in  the  direction  of  the  radar  system. 
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Side  View 
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Figure  5.1:  Grazing  Angle  Error  due  to  Sea  Surface  Slope 

The  effect  of  the  sea  surface  slope  has  a  significant  effect  on  the  local  grazing  angles  for 
high  incidence  radar  look  angles  [refer  to  Figure  5.1  (a)].  However,  for  very  low  incidence 
radar  look  angles,  the  sea  surface  slope,  0r ,  has  small  or  negligible  effects  on  the  local 
grazing  angles,  ^  [refer  to  Figure  5.1  (b)]. 

In  the  composite  sea  clutter  model,  we  did  not  take  into  account  the  correction  to  the 
local  grazing  angle  due  to  sea  surface  slope  when  determining  specular  returns  from  the 
sea  surface  at  near  normal  incidence  angles. 

(a)  Mean-Square  Sea  Slope,  s2 

The  specular  effect  takes  into  account  backscatter  contribution  from  the  large  scale 
waves  or  the  long  gravity  waves.  The  total  mean-square  sea  slope,  s\  is  described 
by  two  components:  (1)  crosswind  (sj!),  and  (2)  up/downwind  (si)]  where  both 
components  are  a  function  of  wind  speed,  U  [2].  An  empirical  model  relating  these 
to  wind  speed  was  determined  from  measurements  obtained  optically  by  Cox  and 
Munk  in  1954  [6]. 

For  a  clean  sea  surface, 

s2c  =  0.003  +  1.92  x  IQ-3U±  0.002  (12) 

s2u  =  0.000  +  3.16  x  10~3f/  ±  0.004 

4  =  s2  +  s2  =  0.003  +  5.12  x  10~3(/  ±  0.004 

For  an  oil  slick  sea  surface, 

s2c  =  0.003  +  0.84  x  10~3U±  0.002  (13) 

si  =  0.005  +  0.78  x  10 ~3U  ±  0.002 

4  =  sc  +  s«  =  0.008  +  1.56  x  10_3(7  ±  0.004 

where  the  ±  value  gives  the  standard  deviation  of  the  observed  value  when  computing 
the  regression  model. 
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Hence  to  determine  the  slope  contribution  of  the  long  gravity  waves  only,  it  is  nec¬ 
essary  to  remove  contributions  to  s2  and  s2  by  the  short  gravity-capillary  waves.  It 
has  also  been  noted  by  Valenzuela  [44]  that  for  sea  backscattering  away  from  normal 
incidence,  only  a  proportion  of  the  total  mean-square  slopes,  Sj.  of  the  sea  surface  is 
included  in  s2  (s2  <  s^),  this  being  the  contribution  of  long  gravity  waves  in  a  fully 
developed  sea  whose  wavelengths  are  greater  than  the  radar  L-band  wavelength. 

Cox  and  Munk  [6]  found  that  the  oil  slick  sea  surface  tends  to  smooth  out  the 
contributions  of  the  capillary  waves.  Therefore,  the  mean  slope  of  the  Cox  and 
Munk  oil  slick  sea  surface  model,  s2  and  s2  of  Equation  (13)  would  be  a  good 
approximation  of  the  proportionate  contribution  of  the  total  mean-square  slope  by 
the  long  gravity  waves.  Valenzuela  [44]  also  obtained  relatively  good  agreement  with 
the  NRL-4FR  measured  radar  backscattered  returns  near  normal  incidence  when  the 
mean-square  slopes  for  a  sea  with  oil  slick  is  used. 

However,  it  has  been  noted  by  Donelan  and  Pierson  [13]  that  using  the  long  gravity 
waves  spectrum  to  determine  the  required  slope  variances,  s2  has  some  shortcom¬ 
ings  in  low  and  high  winds.  The  gravity  waves  on  the  sea  are  usually  higher  than 
the  corresponding  fully  developed  sea  when  winds  are  light  and  lower  during  high 
winds.  Therefore,  instead  of  the  linear  wind  speed  dependence,  we  are  proposing  the 
application  of  a  differential  relationship  as  follows  [6] : 

2  J5.12  x  10_317,  U  <  5m/s  ,  . 

“  1  1.56  x  10-3C/,  U  >  5m/s 


(b)  Reflection  Coefficients,  jR(0) 

Daley  and  Barrick  [44]  found  good  agreement  between  measured  sea  returns  at  near¬ 
normal  incidence  and  predicted  cr^  values  using  Equation  (9).  However,  discrepancies 
between  measured  returns  and  predicted  values  were  observed  for  incidence  angles 
away  from  near-normal.  To  remedy  this,  Valenzuela  suggested  using  the  mean-square 
slopes  derived  for  an  oil  slick  sea  surface  and  adjusting  |i?(0)|  as  follows: 


|R(0) 


0.52 


er  —  1 

(v^+i)2 


(15) 


where  er  is  the  relative  dielectric  constant  of  sea  water  at  L-band  (1.275  GHz)  [13]. 
Note  that  a  constant  0.52  is  used  in  Equation  (15)  instead  of  the  proposed  0.65 
by  Valenzuela  [44].  Better  agreement  between  predicted  backscatter  values,  a°s ,  and 
measured  L-band  JOSH  I  backscatter  values  were  observed  when  using  the  composite 
sea  clutter  model  with  a  constant  of  0.52. 


5.2  Two-Scale  Bragg  Model 

For  grazing  angles  less  than  75°,  the  Bragg  diffraction  model  best  describes  the  diffuse 
backscatter  from  micro- facets  or  short  gravity-capillary  waves  [2].  Guinard  and  Daley  [16] 
proposed  a  composite  sea  surface  scattering  model  based  upon  the  resonant  scattering 
theory  of  the  sea  surface.  The  composite  sea  clutter  model  was  derived  to  remove  the 
constraints  of  the  small  perturbation  method  that  the  product  of  the  vertical  component 
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of  the  wave  number  and  the  surface  displacement  must  be  small.  The  composite  sea 
surface  model  is  based  upon  a  two-scale  model  in  which  a  slightly  rough  surface  (short 
gravity-capillary  waves)  is  superimposed  on  a  larger  wave  structure  (long  gravity  waves). 
The  short  gravity-capillary  wave  with  wave  number,  K  =  2ka  sin  forms  the  facets  for 
radar  backscattering  (also  termed  Bragg  scatterers),  while  the  longer  waves  modulate  the 
Bragg  scatterers  by  “tilting”  or  changing  the  local  incidence  angle  [16]  . 

However,  it  has  been  observed  by  Valenzuela  and  associates  [46]  that,  the  backscattered 
sea  clutter  for  vertical  polarisation  is  insensitive  to  the  tilting  effect.  Therefore  for  practical 
purposes,  the  backscatter  signal  for  vertical  polarisation  sea  clutter  can  be  predicted  by 
slightly  rough  surface  scattering  theory  with  no  modification  to  the  radar  depression  angle 
except  to  account  for  diffraction  due  to  non- homogeneous  atmospheric  medium. 


5.2.1  First-Order  Composite  Sea  Clutter  Model  (CSM),  v°Bi(pp) 

The  composite  sea  clutter  model,  hereinafter  referred  to  as  the  CSM  sea  clutter  model 
is  a  two-part  first-order  model,  whereby  the  scattering  from  any  composite  rough  sea 
surface  whose  wave  length  is  Bragg  resonant  and  height  is  given  by  Phillips’  proposed 
equilibrium  range  spectrum  [16]. 

Horizontal  Polarisation,  <?°b\(Hh)  —  4-7T  k4  sin4  ip  olBh  W(Kx,Ky)  (16) 

■v"  . 

Surface  Height  Spectrum 

Vertical  Polarisation,  a°Bi  =  Air  k4  sin4  ip  ayy  W(Kx,Ky) 

Surface  Height  Spectrum 

where  k0  is  the  radar  wave  number,  ip  is  the  grazing  angle,  and  Kx,  Ky  are  the  wave 
numbers  in  the  x  and  y  directions  on  the  infinite  surface  [16]. 

Philips  found  that  at  equilibrium  range,  the  L-band  sea  spectrum  has  the  form 

W{K)  »  BK~4  (17) 


where  K  =  +  K$  and  B  =  6  x  10  3  for  gravity  waves. 

The  CSM  model  for  an  L-band  radar  system  could  therefore  be  simplified  to  the  fol¬ 
lowing: 

Horizontal  Polarisation,  o°Bi^HH)  =  1-5  x  10~37t  olhh  tan4  ip  (18) 

Vertical  Polarisation,  cr Bi(w)  =  1-5  x  10_37t  ctvv  tan4  ip 


where  ip  is  the  local  grazing  angle. 

The  first-order  scattering  coefficient,  app  is  [16] 


Horizontal  Polarisation,  olhh  =  - .  .  ^ -  - 

[cos  0i  +  y  (er  -  shV  6i)Y 

„  ,.  ,  r,  ,  ■  ..  (er  —  l)[er(l  +  sin2  OP)  —  sin2  0,] 

Vertical  Polarisation,  ayv  =  - ,  „  =— — 

[ercos6>t  +  v/^^irt)]2 


(19) 
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Airborne  Sea  Clutter  for  L-band  VV  Polarisation 


where  er  is  the  relative  dielectric  constant  of  the  sea  surface  with  value  72  —  j 59  at  L-band 
(1.275  GHz)  [13], 

A  comparison  of  the  CSM  sea  clutter  model  with  the  NRL-4FR  data  recorded  in  the 
upwind  direction  with  sea  surface  wind  speed  of  24  m/s  and  swell  of  8  meters  is  shown  in 
Figure  5.2: 

1.  The  CSM  model  provides  a  realistic  upper  bound  for  VV-polarisation  data  for  grazing 
angle,  ip  <  75°; 

2.  To  the  first  order,  cr°Bl^vv^  is  given  by  the  spectral  density  of  the  resonant  Bragg 
waves  [51]. 

The  CSM  sea  clutter  model  estimates  the  backscattering  from  the  sea  surface  under 
saturation  conditions.  The  model  assumes  the  spectral  density  of  the  resonant  Bragg 
wave  to  be  independent  of  wind  speed  under  saturation  conditions  [33].  As  shown  in 
Figure  5.2,  the  CSM  sea  clutter  model  provides  a  realistic  sea  clutter  model  for  high 
sea  state  conditions,  e.g.,  for  wind  speed  of  23.5  m/s  which  according  to  Nathanson  [29] 
is  categorised  as  strong  gale  condition.  However,  a  more  appropriate  sea  clutter  model 
responsive  to  wind  speed  changes  is  required  to  evaluate  the  sea  backscattering  coefficient, 
<r°,  for  various  sea  conditions. 
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5.2.2  Second-Order  Composite  Sea  Clutter  Model  ( MCSM  ) 


The  original  CSM  model  was  modified  into  a  three-part,  second-order  sea  clutter  model 
and  hereafter  referred  to  as  the  Modified  Composite  Sea  Clutter  model  or  MCSM  in  short. 
The  MCSM  sea  clutter  model  is  based  upon  the  concept  of  slightly  rough  sea  surface 
scattering  whose  length  is  Bragg  resonant  and  height  is  given  by  a  second-order  directional 
sea  spectrum  proposed  by  Fung  [14]  and  Chen  [5].  The  model  also  accounts  for  specular 
return  near-normal  incidence 

MCSM  Modified  Sea  Clutter  Model 

Reflectivity  equation  (dB  m2/m2), 

a°  =  a°S  +  a°B2(pp)  (20) 


1.  Second  Order  Bragg  Scattering  Model,  o°B 2 ^ 

Horizontal  Polarisation,  a°B2{HH)  =  k^sm4ipaHH  W (K,  <j) )  (21) 

Directional  Sea  Spectrum 

Vertical  Polarisation,  =  fc2  sin4  ^  av'V'  W(K,<j>) 

Directional  Sea  Spectrum 

where  k0  is  the  radar  wavenumber,  ip  is  the  grazing  angle,  K  is  the  sea  wavenumber, 
and  4>  is  the  wind  direction  with  respect  to  boresight. 


2.  The  first-order  scattering  coefficient,  app  is  [16] 
Horizontal  Polarisation,  olhh  = 


(er  -  1) 


Vertical  Polarisation,  ayy 


[cos  9i  +  \J (f-T  -  sin2  d,)]2 
(er  -  l)[er(l  +  sin2^)  -  sin2  0j] 
[er  cos  9i  -I-  i / (er 


(22) 


sin2  9i)]2 

where  er  is  the  relative  dielectric  constant  of  the  sea  surface  with  value  72  —  j59 
(L-band,  1.275  GHz)  [13]. 


3.  Two-dimensional  (2-D)  Sea  Ripple  Spectrum 

Adopting  the  2-D  sea  ripple  spectrum  developed  by  Pierson  [43]  and  modified  by 
Fung  and  Lee  [5],  the  2-D  directional  (taking  into  account  the  wind  direction)  sea 
spectrum  is  [43,  5] 

W{K,(t>)  =  -!-(l  +rcos2  4)W(K)  (23) 


(a)  The  constant,  r,  is  given  by 


where  v  is  the  ratio  of  the  slope  variances  in  the  crosswind  and  upwind  directions 
[15,  43,  5], 


0.003+1. 92xl0~3t/ 
3.16x10 -*U  ’ 

0.003+0. 84xl0~3t/ 
0. 005+0. 78xlO-3U’ 


U  <  5m/s 
U  >  5m/s 


and  U  (m/s)  is  the  wind  speed  at  an  altitude  of  12.5  m  above  the  sea  surface. 
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(b)  Wind  Direction,  4> 


Wind  direction,  </> 


'0°, 
90°or 
k 180°, 


90° 


upwind 

crosswind 

downwind 


(c)  1-D  sea  ripple  spectra  for  a  single  wave  train,  W(K)  [15,  43,  5] 

(i-p) 

(1+3#^). <5 

W{K)  =  A(  2tt)p~1 


9 


ro.922  ^  +  Jj)' 


(1+P) 

2 


(24) 


where  A  =  1.75  if  U  <  5  m/s  or  A  =  0.875  if  U  >  5  m/s,  K ^  =  ^  = 
(3.63)2cm~2,  g  =  acceleration  of  gravity  9.81  m/s2,  p  =  sea  water  density,  and 
r  =  sea  surface  tension.  Note  that  according  to  Valenzuela  [45],  K  should  have 
a  —3.922  power  law  rather  than  the  usual  —4  law  for  L-band  resonant  Bragg 
waves.  This  has  been  adopted  in  the  MCSM  model. 

According  to  Fung  [15]  and  Chen  [5]  the  relationship  between  the  sea  surface 
wind  speed  at  an  altitude  z  and  friction  velocity  at  neutral  stability  of  the 
atmosphere,  U*  is  [17] 

u.&mU-) 

0.4  \Z0J 

where 


and  p  =  5 


Zq  (cm)  = 

logio  U*  ■ 


0.684  \ 
U*  ) 


+  4.28  x  10 ~5U*2  -  0.0443 


5.2.3  The  Use  of  the  NRL-4FR  Data  to  Test  the  CSM  and  MCSM 
Models 

To  check  the  validity  of  the  models,  we  compared  the  CSM  and  MCSM  backscatter 
predictions  with  the  NRL-4FR  data  [46].  Plots  of  the  predicted  and  measured  upwind 
backscattering  coefficients,  a°,  versus  grazing  angles  for  varying  sea  surface  wind  conditions 
are  shown  in  Figures  5.3  (a)  -  (e). 

The  plots  of  Figures  5.3  (a)  -  (e),  show  a  slight  increase  in  <r°  with  an  increase  in 
wind  speed.  Good  agreement  between  predicted  backscatter  values  and  measured  data  is 
observed  for  intermediate  upwind  wind  speed,  5  <  wind  speed  <  20  m/s  for  the  MCSM  sea 
clutter  model.  However,  the  CSM  model  overestimates  the  backscattering  coefficients,  cr°, 
for  upwind  speed  within  this  intermediate  range.  A  similar  result  is  observed  for  very  low 
wind  speed  (2.5  m/s)  but  good  agreement  between  predicted  CSM  values  and  measured 
data  is  observed  for  very  high  wind  speed  (23  m/s). 

At  very  low  wind  speed  (2.5  m/s),  the  MCSM  model  also  over-estimates  the  upwind 
backscattering  coefficient,  cr°,  for  low  grazing  angles,  but  good  agreement  is  observed  at 
higher  grazing  angles.  However,  the  opposite  is  observed  at  very  high  surface  wind  speed 
condition  (23.5  m/s).  Disagreement  between  the  predicted  MCSM  and  measured  values 
becomes  evident  at  higher  grazing  angles  for  very  high  wind  speed  (23.5  m/s). 
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6  Comparisons  with  Radar  Measurements 

As  discussed  earlier  in  Section  4.1,  poor  agreement  between  the  TSC  and  HYB  sea 
clutter  models  and  the  measured  NRL-4FR  data  were  observed  for  two  wind  speeds  (5.6 
m/s  and  12.1  m/s).  We  revisited  the  exercise  of  testing  the  accuracy  of  these  semi- 
empirical  models  with  the  same  measured  NRL-4FR  data,  but  this  time  we  included  the 
MCSM  model  in  the  exercise. 


Figure  6.1:  Dependence  of  G°m  on  Wind  Speed ,  Vw  ( L-band ,  VV-Pol.  and  Upwind),  (a) 
Wind  speed  =  5.6  m/s;  (b)  Wind  speed  =  12.1  m/s; 


Plots  of  the  results  in  Figures  6.1  (a)  and  (b)  show  good  agreement  between  the 
measured  NRL-4FR  data  and  the  MCSM  predicted  a°  values  for  ip  >  3°  for  the  two  wind 
speeds.  For  ip  <  3°,  good  agreement  between  Nathanson’s  averaged  data  and  MCSM 
predicted  o°  values  occurs  at  low  wind  speed  (5.6  m/s),  but  discrepancies  occur  at  higher 
wind  speed  (12.1  m/s). 

We  further  tested  the  validity  of  the  MCSM  sea  clutter  model  in  predicting  L-band 
V-pol.  sea  clutter  under  varying  sea  conditions  using  the  JOSH  I,  JOSS  II,  and  the 
POLSAR  DARWIN  measured  data.  Figures  6.2  -  6.4  (a)  -  (d)  again  show  the  poor 
agreement  between  the  TSC  and  HYB  models  with  the  JOSS  I  measured  data  for  varying 
wind  speeds  and  radar  look  directions.  Where  these  two  models  fail,  the  MCSM  shows 
good  agreement  with  the  measured  data  for  all  wind  speeds,  radar  look  directions  and 
grazing  angles.  Similar  results  were  observed  with  the  JOSS  II  measured  data,  shown  in 
Figures  6.5  -  6.7  (a)  -  (c). 

For  the  POLSAR  data  acquired  over  the  DARWIN  region  under  low  wind  speed  con¬ 
dition,  the  MCSM  model  again  accurately  predicted  the  backscattering  coefficients,  cr°, 
for  20°  <  ip  <  65°.  Figure  6.8  shows  the  measured  data,  and  the  predicted  a°  values  by 
the  MCSM  and  TSC  sea  clutter  models.  The  TSC  model  underestimates  the  backscatter 
values  when  compared  to  the  measured  data. 

We  have  tested  four  sea  clutter  models  against  measured  data  obtained  from  calibrated 
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Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


JOSS  I  (Wind  speed  =  5  8  m/s) 
•  GIT  (Wind  speed  =  5  8  m/s) 
HYB 
TSC 

■  MCSM  (Wind  speed  =  5  8  m/s) 


*  JOSS  I  (Wind  speed  =  84  m/s) 
-  -  GIT  (Wind  speed  =  8  4  m/s) 

-  HYB 

-  MCSM  (Wind  speed  «  8  4  rVs) 


Grazing  Angle,  \j /  (degrees) 


Grazing  Angle,  v  (degrees) 


Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


JOSS  I  (Wind  speed  =  1 1 .55  m/s) 
-  GIT  (Wind  speed  =  1 1.55  m/s) 
HYB 

•  MCSM  (Wind  speed  =  1 1  55  m/s) 


Grazing  Angle,  v  (degrees) 


Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


•  JOSS  I  (Wind  speed  =  13  125  m/s) 

- GIT  (Wind  speed  =  13.125  m/s) 

-  HYB 

-  MCSM  (Wind  speed  =  13  125  nVs) 


Grazing  Angle,  i|/  (degrees) 


Figure  6.2:  Dependence  of  o°  on  wind  speed  -  JOSH  I,  U  (L-band,  VV-Pol.  and  upwind), 
(a)  wind  speed  =  5.8  m/s;  (b)  wind  speed  =  8-4  m/s;  (c)  wind  speed  =  11.55  m/s;  (d) 
wind  speed  =  13.125  m/s 
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Figure  6.3:  Dependence  of  a°  on  wind  speed  -  JOSH  I,  U  (L-band,  VV-Pol.  and  down¬ 
wind).  (a)  wind  speed  =  5.8  m/s;  (b)  wind  speed  =  8.4  m/s;  (c)  wind  speed  =  11.55  m/s; 
(d)  wind  speed  =  13.125  m/s 
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Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


(b) 


Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


Airborne  Sea  Clutter  for  L-band  (VV-Pol.) 


(d) 


Figure  6-4-'  Dependence  of  a0  on  wind  speed  -  JOSH  I,  U  (L-band,  VV-Pol.  and  cross- 
wind).  (a)  wind  speed  =  5.8  m/s;  (b)  wind  speed  —  8.4  m/s;  (c)  wind  speed  =  11.55  m/s; 
(d)  wind  speed  —  13.125  m/s 
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Airborne  Sea  Clutter  for  L-band  (VV-Pol) 


systems  under  varying  sea  conditions,  radar  look  directions  and  grazing  angles.  By  cover¬ 
ing  a  wide  range  of  sea  conditions,  we  were  able  to  select  the  most  appropriate  sea  clutter 
model.  These  results  suggested  that  among  these  four  ( GIT  ,  HYB  ,  TSC  ,  MCSM )  sea 
clutter  models,  the  MCSM  model  is  the  better  generic  sea  clutter  model  for  predicting  the 
sea  surface  backscattering  coefficient,  a°,  for  an  L-band,  VV-polarised  radar  system. 

We  have  not  considered  the  inclusion  of  contributions  from  wedge  scattering  proposed 
by  Lyzenga  and  associates  [24],  due  to  several  reasons.  Firstly,  it  is  impossible  to  mea¬ 
sure  the  number  of  wedges  per  square  meter,  and  their  average  length.  Secondly,  it  has 
been  noted  by  Lyzenga  that  for  vertical  polarisation,  there  is  good  agreement  between 
Bragg  scattering  and  the  NRL-4FR  data  with  minimal  contribution  from  wedge  scatter¬ 
ed.  Thirdly,  at  L-band,  the  wedge  appear  to  reach  saturation  and  break  slower  compare  to 
X-band  wedges.  In  view  of  this,  the  contributions  from  Bragg  scatterers  can  be  extended 
to  low  grazing  angles  for  L-band  VV-polarised  system. 
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7  Radar  Backscattering  Signals  from  Darwin 

Coastal  Areas 

For  a  radar  system  operating  close  to  the  coastal  region,  its  detection  and  tracking  ca¬ 
pabilities  will  be  severely  affected  by  compounded  backscatter  reflections  from  the  coastal 
water,  foreshore  vegetation,  and  man-made  structures.  The  clutter  return  from  the  coastal 
land  areas  will  depend  on  the  coastal  terrain  characteristics  e.g.,  vegetation  type,  presence 
of  man-made  structure,  ground  surface  roughness  and  many  more. 

Calibrated  POLSAR  images  of  the  Darwin  areas,  CM5512  and  CM5148  shown  in 
Figures  3.4  and  3.5  were  used  to  study  the  properties  of  mean  backscattering  coefficient, 
cr°,  from  coastal  water  and  terrain.  The  availability  of  the  Darwin  POLSAR  data  allows 
us  to  study  the  characteristics  of  backscatter  signals  from  coastal  water,  coastal  saline  flat 
and  mangrove  forest  and  their  dependence  on  grazing  angle  for  an  L-band  VV-polarised 
radar  system.  Terrain  scatterer  types  in  the  POLSAR  images  were  identified  using  the 
TOPO-IOOK  and  TOPO-50K  maps. 


7.1  Coastal  Water  Backscattering  Characteristics 

The  radar  signal  backscatter  from  the  coastal  water  surface  close  to  the  coastal  saline 
flat  of  the  Darwin  coast  at  low-tide  were  measured  from  the  POLSAR  images.  As  shown 


Airborne  Shallow  Water  Clutter  for  L-band  (VV-Pol) 


Figure  7.1:  Radar  Backscattering  Coefficients  from  Shallow  Water  Surface 

in  Figure  7.1,  the  MCSM  sea  clutter  model  has  good  agreement  with  the  measured  coastal 
water  backscattered  values.  However,  with  increasing  grazing  angles  the  returns  increases 
more  rapidly  than  that  predicted  by  the  MCSM  sea  clutter  model. 
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7.2  Wetland  Features 


Radar  backscatter  from  wetland  is  dependent  on  the  tide  conditions.  At  high  tide, 
the  sea  water  will  completely  cover  the  aerial  roots,  resulting  in  a  reduction  in  multipath 
backscatter  aerial-ground  interaction  (D)  shown  in  Figure  7.2.  At  the  same  time,  flooded 
mangrove  forest  have  a  smoother  surface  than  a  non-flooded  surface.  Radar  backscatter 
from  a  flooded  surface  tends  to  be  higher  by  3  -  10  dB  compared  to  non-flooded  surface 
[28].  L-band  signals  are  able  to  penetrate  the  vegetation  canopy  to  cause  either  direct  or 
multipath  reflection  on  the  smooth  water  surface  (A,C,D). 


Canopy  Layer 


Trunk  Layer 
.  Aerial  root  layer 


Figure  7.2:  Backscatter  components  of  wetlands.  A  =  surface  backscatter.  B  =  canopy 
backscatter.  C  =  double  bounce  trunk-surface  backscatter.  D  =  double  bounce  aerial- 
surface  backscatter.  E  =  double  bounce  canopy- surface  backscatter.  [48] 


7.2.1  Radar  Backscattering  from  Coastal  Saline  Flat  Surfaces  (A) 

L-VV  backscattered  values  from  coastal  saline  flat  surfaces  under  non-flooded  condi¬ 
tions  were  evaluated  from  the  POLSAR  images  of  Figures  3.4  and  3.5  for  varying  graz¬ 
ing  angles.  The  radar  backscattering  characteristic,  cr°,  dependence  on  grazing  angles 
for  coastal  saline  flat  is  shown  in  Figure  7.3.  The  result  shows  an  increasing  trend  in 
backscatter  with  increasing  grazing  angles. 

The  Small  Perturbation  Model  (SPM)  has  been  used  to  calculate  the  backscattering 
and  specular  reflection  from  the  coastal  saline  flat  surfaces  for  an  L-band  radar.  The 
backscattering  coefficient,  cr°pp,  is  given  by  [43,  11] 

<o  =  8 {khrmsf  cos4  6i.  W{2kosm0i,0\  \app\2  (25) 

^  V 

Roughness  Spectrum 


where 

k0 

=  radar  wavenumber  (k0  =  x); 

h-rms 

=  surface  rms  roughness  (m); 

L 

=  surface  correlation  length  (m 

Oi 

=  incidence  angle  (radian); 

&pp 

=  scattering  coefficient. 

Note  that  the  SPM  method  is  valid  only  if  khrms  <  0.3  and  kL  <  3  [31]. 
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Darwin  Coastal  Saline  Flat  Reflectivity  Model  (L-band,  VV) 


Grazing  Angle  (degrees) 

Figure  7.3:  Radar  Backscattering  from  Costal  Saline  Flat 


1.  The  first  order  scattering  coefficient,  app  [16] 

,  „  ,  .  .  (Cr  ~  1) 

Horizontal  Polarisation,  olhh  =  - T,  . 

[cos  9i  +  y/(er  -sm20i)]2 

(er  -  l)[er(l  +  sin2  Of)  -  sin2  Oj 


Vertical  Polarisation,  avv  =  r  „  r, - .  o  ^ 

[er  cos  di  +  y  (er  -  sinz  Oi)Y 

where  er  is  the  relative  dielectric  constant  of  the  coastal  saline  flat  surface. 


2.  Surface  Roughness  Spectrum 

Gaussian  Correlation  Coefficient  [31,  43] 

For  a  lsi  order  Gaussian  correlation  coefficient  of  the  form,  p(£)  =  exp(-£2/ L2),  the 
isotropic  surface  roughness  spectrum  is  given  by  the  zero-order  Hankel  transform, 

W(2kosmdi,0)  =  Ho{p(£)} 

=  exp[—  (kL  sin  di)2] 


hrms  =  0.033  cm  and  L  =  9.2  cm,  and  a  dielectric  constant  of  72  -  j59  for  the  coastal 
saline  fiat  surface  were  assumed  in  the  simulation.  The  basis  for  these  assumptions  lies  in 
the  reality  of  the  coastal  saline  flat  surface  conditions:  (1)  the  area  is  very  level,  and  (2) 
because  of  the  twice  daily  inundation  of  sea  water  at  high-tide,  the  non-flooded  coastal 
saline  flat  surface  is  saturated  with  sea  water.  Therefore  the  dielectric  property  of  the 
coastal  saline  flat  surface  is  similar  to  that  of  the  sea  water  [48].  However,  the  dielectric 
content  will  vary  during  the  monsoon  rainy  season  due  to  the  dilution  of  the  salt  water. 
Figure  7.4  shows  the  measured  coastal  saline  flat  data,  and  the  predicted  <jSvv  values  by 
the  SPM  clutter  model  with  Gaussian  correlation  coefficient. 
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Airborne  Coastal  Saline  Flat  Clutter  for  L-band  (VV) 


Figure  7.4'-  L-VV  Backs cattering  Prediction  for  Non-flooded  Coastal  Saline  Flat  Surfaces 


Comparison  of  the  predicted  backscattering  with  the  POLSAR  data  shows  good  agree¬ 
ment  for  returns  from  non-flooded  coastal  saline  flat  surfaces.  However,  due  to  the  limited 
data  available,  further  validation  would  have  to  be  carried  out  to  assess  the  suitability  of 
the  SPM  model  in  predicting  backscatter  from  non-flooded  coastal  saline  flat  surface. 


7.2.2  Radar  Backscattering  from  Coastal  Mangrove  Forest 

The  double-bounce  ground  trees/  aerial  roots  backscatter  signal  is  enhanced  by 
the  presence  of  sea  water  under  the  Rhizophora  stylosa  trees.  The  strong  double-bounce 
backscattered  signal  is  observed  along  the  unshadowed  foreshore  section  of  the  coastal 
Darwin  areas.  Mouchot  and  Garello  [27]  also  noted  that  the  presence  of  vertical  structure 
of  vegetation  in  both  swamps  and  mangrove  forest  relative  to  the  surface  of  the  sea  acts 
as  a  corner  reflector,  resulting  in  an  area  of  bright  return  for  L-band  radars. 

Strong  returns  can  be  seen  from  mangrove  forest  surrounding  the  coastal  Darwin  areas 
in  Figures  3.4  and  3.5.  Richard  and  associates  [36]  have  shown  that  the  specular  returns 
from  tree  trunks  to  ground  and  then  back  to  the  radar  system  give  rise  to  higher  scattering 
levels  in  comparison  to  volume  scattering  from  the  canopy  or  returns  from  the  ground 
surface.  Measured  backscatter  signal  from  the  mangrove  forest  surrounding  the  coastal 
Darwin  areas  is  shown  in  Figure  7.5  to  increase  with  grazing  angles.  This  could  be  partly 
explained  by  the  shorter  path  length  through  the  canopy  and  trunk  layers  at  higher  grazing 
angles. 

Models  relating  the  structure  and  geometry  of  the  mangrove  forest  were  not  carried 
out  due  to  the  lack  ground  truth  information. 
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Darwin  Mangrove  Forest  Reflectivity  Model  (L-band,  VV) 


Grazing  Angle  (degrees) 


Figure  7.5:  Radar  Backscattering  from  Mangrove  Forest  (L-VV) 

8  Conclusions 

Data  obtained  from  the  NRL  4FR  system  (NRL-4FR,  JOSS  I,  JOSS  II),  Nathanson 
compiled  measurements,  and  the  POLSAR  Darwin  datasets  have  been  collated  with  five 
sea  clutter  models  (GIT  ,  HYB  ,  TSC  ,  GSM  ,  and  MCSM ).  These  datasets  were  used  to 
determine  the  validity  of  the  five  sea  clutter  models  in  predicting  the  backscatter  signals 
from  the  sea  surface  for  varying  sea  conditions,  radar  look  directions  and  grazing  angles. 

The  result  of  the  comparisons  have  shown  that  the  two-scale  model  modified  for 
anisotropic  sea  spectrum,  MCSM  ,  accurately  predicted  the  NRL-4FR,  JOSS  I,  and  JOSS 
II  data  for  a  wide  range  of  grazing  angles,  look  angles,  and  sea  conditions.  In  addition,  the 
MCSM  sea  clutter  model  accurately  predicted  the  measured  backscattering  coefficients, 
cr°,  for  sea  surface  around  the  Port  Darwin  area  for  grazing  angles  between  25°  to  65°. 
However,  this  model  will  have  to  be  verified  with  backscattered  signal  acquired  near  normal 
incidence  from  sea  temperature  around  30°  C. 

Note  that  the  work  in  this  report  has  been  primarily  aimed  at  L-band  VV-polarised 
radar  system.  Further  research  will  have  to  be  carried  out  to  assess  the  suitability  of 
applying  the  MCSM  model  to  predict  the  sea  backscattering  coefficient,  o  ,  for  varying 
sea  conditions  for  radar  systems  operating  in  other  frequency  bands. 
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Appendix  A:  POLSAR  Mean  Backscattering 
Coefficients,  a°  (Darwin  Data) 


The  data  in  this  section  are  organised  by  sea  and  terrain  categories  in  Table  A1  to  A4. 
For  each  category,  the  table  provides  a  statistical  summary  of  the  measured  L-band  VV- 
polarised  backscattered  POLSAR  Darwin  data.  At  each  grazing  angle  for  which  data  are 
measured,  the  following  information  is  provided  in  the  statistical  distribution  tables. 

1.  ip,  the  grazing  angle  in  degrees 

2.  N,  the  number  of  independent  samples  of  pixels  used  in  the  measurement  process 
that  led  to  each  cr°. 

3.  cr°,  the  normalised  mean  backscattering  coefficients  in  dB. 

4.  The  median,  standard  deviation,  minimum,  and  maximum  values  of  the  a°  distri¬ 
bution,  expressed  in  dB. 
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Table  Al:  Backscattering 

Coefficient  for 
Median  (dB) 

Deep  Water  for  L- 

VV  POLSAR  Darwin 

N  Samples 

a°  (dB) 

Std.  Dev.  (dB) 

Min  (dB) 

Max  (dB) 

28 

828 

-24.96 

-25.50 

2.55 

-34.50 

-19.00 

29 

54648 

-22.69 

-23.00 

2.06 

-35.25 

-12.50 

30 

56942 

-22.07 

-22.50 

2.24 

-37.50 

-15.00 

31 

76018 

-22.18 

-22.50 

2.34 

-39.25 

-15.50 

32 

77080 

-22.40 

-22.75 

2.33 

-40.00 

-15.75 

33 

63390 

-22.23 

-22.50 

2.84 

-40.00 

-15.25 

34 

41400 

-21.40 

-21.75 

2.06 

-39.50 

-15.00 

35 

38916 

-21.39 

-21.75 

1.99 

-33.00 

-6.75 

36 

45284 

-21.17 

-21.50 

2.13 

-34.50 

-5.75 

37 

43900 

-21.09 

-21.50 

2.20 

-34.00 

-5.50 

38 

33948 

-20.01 

-21.25 

2.70 

-36.00 

7.25 

39 

32292 

-18.21 

-21.00 

2.96 

-33.25 

11.75 

40 

31050 

-19.83 

-20.50 

2.41 

-35.00 

3.25 

41 

23184 

-20.15 

-20.50 

2.26 

-33.50 

-11.75 

42 

7644 

-19.63 

-20.00 

2.30 

-35.25 

-13.75 

43 

11368 

-19.47 

-19.75 

2.37 

-36.25 

-12.00 

44 

10976 

-19.12 

-19.50 

2.47 

-36.50 

-2.25 

45 

10388 

-19.66 

-20.00 

3.00 

-36.50 

-12.75 

48 

3838 

-18.89 

-19.25 

1.91 

-27.50 

-13.75 

49 

4747 

-17.73 

-18.00 

1.87 

-25.50 

-10.00 

50 

4646 

-17.44 

-17.75 

1.89 

-26.50 

-11.50 

51 

8532 

-17.84 

-18.25 

1.90 

-27.75 

-12.25 

52 

10722 

-17.92 

-18.25 

2.52 

-35.25 

-7.50 

53 

22972 

-17.54 

-18.25 

3.85 

-40.00 

-8.75 

54 

24477 

-16.57 

-17.00 

3.70 

-40.00 

0.00 

55 

35139 

-16.77 

-17.25 

3.09 

-40.00 

0.00 

56 

39671 

-16.28 

-16.75 

2.88 

-40.00 

0.00 

57 

32984 

-16.21 

-16.75 

2.67 

-40.00 

0.00 

58 

36026 

-15.82 

-16.25 

2.46 

-36.75 

-8.75 

59 

46881 

-15.51 

-16.00 

2.63 

-35.25 

-7.00 

60 

48099 

-15.44 

-15.75 

2.93 

-36.75 

-4.25 

61 

47190 

-15.16 

-15.50 

2.90 

-36.00 

-3.50 

62 

47020 

-14.90 

-15.25 

2.76 

-35.50 

-2.00 

63 

45906 

-14.16 

-14.50 

2.47 

-34.50 

-3.50 

64 

44698 

-13.67 

-14.00 

2.46 

-30.50 

-6.75 

65 

43584 

-12.92 

-13.25 

2.58 

-28.00 

-6.00 

66 

44206 

-11.85 

-12.25 

2.20 

-40.00 

-5.50 

67 

21211 

-10.90 

-11.25 

2.36 

-40.00 

-4.25 
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Table  A2:  Backs  cattering  Coefficient  for  Coastal  Water  for  L-VV  POLSAR  Darwin 


i> 

N  Samples 

a°  (dB) 

Median  (dB) 

Std.  Dev.  (dB) 

Min  (dB) 

Max  (dB) 

29 

4087 

-21.69 

-22.00 

1.87 

-29.25 

-16.75 

30 

8701 

-21.70 

-22.00 

2.64 

-40.00 

-15.75 

31 

17413 

-21.79 

-22.25 

2.82 

-39.25 

-14.50 

32 

6250 

-20.89 

-21.25 

2.52 

-37.25 

-14.00 

33 

3625 

-21.61 

-22.00 

2.84 

-38.00 

-15.25 

34 

2418 

-21.57 

-22.25 

2.48 

-33.75 

-15.25 

35 

10130 

-20.08 

-20.75 

2.82 

-35.25 

-13.75 

36 

22201 

-20.34 

-20.75 

2.38 

-34.00 

-12.00 

37 

11708 

-21.32 

-21.75 

2.45 

-39.25 

-14.50 

38 

6042 

-21.47 

-21.75 

2.32 

-35.25 

-15.25 

39 

11346 

-20.57 

-21.00 

2.22 

-33.50 

-14.00 

40 

7093 

-20.31 

-20.75 

2.43 

-32.75 

-14.25 

41 

4533 

-19.87 

-20.25 

2.14 

-32.50 

-13.50 

42 

6473 

-20.11 

-20.50 

2.13 

-32.00 

-15.00 

43 

3087 

-19.38 

-19.75 

2.03 

-29.75 

-13.00 

44 

8079 

-19.33 

-19.75 

2.04 

-29.00 

-13.00 

45 

7534 

-19.11 

-19.50 

1.88 

-28.75 

-14.00 

46 

5784 

-19.14 

-19.50 

1.98 

-29.75 

-12.75 

47 

4544 

-19.48 

-19.75 

1.94 

-29.75 

-14.75 

48 

1533 

-19.23 

-19.50 

1.92 

-26.50 

-12.00 

49 

3068 

-18.53 

-18.75 

1.88 

-25.75 

-14.25 

50 

4692 

-18.46 

-18.75 

1.96 

-28.75 

-12.25 

51 

3960 

-17.60 

-18.00 

1.90 

-25.75 

-11.50 

52 

3990 

-17.50 

-17.75 

2.10 

-26.50 

-10.25 

53 

3876 

-16.87 

-17.25 

1.98 

-26.50 

-11.25 

54 

3923 

-17.32 

-17.75 

2.06 

-25.00 

-11.75 

55 

4606 

-17.68 

-18.00 

1.92 

-25.50 

-10.75 

56 

1549 

-17.76 

-18.00 

1.84 

-24.75 

-13.75 

57 

3555 

-16.82 

-17.00 

1.80 

-26.00 

-11.00 

58 

3397 

-15.45 

-15.75 

1.84 

-24.75 

-10.75 

59 

2636 

-14.22 

-14.50 

1.96 

-22.00 

-8.50 

60 

6399 

-13.63 

-14.00 

2.01 

-23.50 

-7.00 

61 

5661 

-13.88 

-14.25 

1.89 

-23.75 

-8.00 

62 

3567 

-13.97 

-14.25 

1.82 

-21.00 

-8.75 

63 

5418 

-13.70 

-14.00 

1.85 

-22.25 

-7.75 

64 

7536 

-13.03 

-13.25 

1.85 

-25.50 

-8.00 

65 

5437 

-11.83 

-12.50 

2.20 

-20.00 

-5.50 

66 

930 

-9.72 

-10.00 

1.81 

-17.50 

-5.50 

67 

840 

-9.43 

-9.75 

1.80 

-16.75 

-5.00 

43 
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Table  A3:  Backs cattering  Coefficient  for  Coastal  Saline  Flat  for  L-  VV  POLSAR  Darwin 


N  Samples 

<j°  (dB)  : 

Median  (dB) 

Std.  Dev.  (dB) 

i  Min  (dB) 

Max  (dl 

26 

11178 

-23.49 

-26.75 

4.81 

-40.00 

-8.50 

27 

19238 

-25.11 

-26.00 

2.77 

-36.75 

-14.75 

28 

1859 

-25.35 

-29.00 

5.90 

-40.00 

-15.50 

29 

28523 

-23.42 

-24.75 

4.24 

-40.00 

-9.75 

30 

15322 

-23.07 

-24.75 

4.11 

-38.75 

-7.50 

31 

21869 

-23.30 

-24.50 

3.37 

-36.25 

-13.25 

32 

40784 

-25.70 

-27.75 

3.53 

-38.00 

-12.50 

33 

3646 

-22.82 

-24.00 

3.77 

-36.50 

-14.00 

34 

5404 

-20.52 

-21.25 

2.54 

-31.25 

-13.50 

35 

4808 

-19.56 

-21.00 

2.85 

-31.25 

-11.00 

36 

9121 

-23.40 

-24.25 

3.53 

-37.50 

-12.25 

37 

10898 

-24.32 

-25.25 

3.35 

-36.25 

-17.75 

38 

10217 

-23.69 

-24.25 

3.39 

-37.50 

-16.00 

39 

506 

-22.01 

-22.75 

2.52 

-30.25 

-14.50 

40 

744 

-20.29 

-21.25 

2.77 

-29.25 

-13.50 

42 

1722 

-19.91 

-20.75 

2.57 

-29.50 

-11.25 

43 

2860 

-19.89 

-20.50 

2.47 

-31.25 

-13.25 

44 

1707 

-18.96 

-19.75 

2.55 

-28.25 

-10.50 

45 

486 

-19.16 

-20.00 

2.59 

-29.50 

-13.00 

46 

3851 

-21.55 

-22.50 

2.69 

-32.25 

-12.75 

47 

9378 

-18.03 

-19.00 

3.58 

-32.75 

-9.75 

48 

10830 

-16.73 

-18.25 

3.87 

-37.75 

-6.75 

49 

15697 

-16.13 

-18.75 

4.72 

-36.50 

-4.50 

50 

15580 

-15.91 

-18.00 

4.57 

-37.75 

-3.50 

51 

18824 

-15.17 

-17.00 

4.10 

-34.50 

-4.50 

52 

13614 

-14.61 

-16.25 

4.18 

-36.00 

-3.75 

53 

17721 

-15.65 

-17.50 

4.21 

-37.25 

-3.50 

54 

16459 

-15.50 

-17.00 

4.62 

-39.50 

-5.00 

55 

13096 

-17.52 

-19.00 

3.25 

-31.00 

-5.75 

56 

12040 

-18.04 

-20.25 

3.57 

-29.75 

-5.00 

57 

10179 

-17.43 

-19.25 

3.31 

-29.75 

-4.75 

58 

4827 

-14.84 

-15.50 

2.99 

-30.25 

-7.00 

59 

196 

-13.00 

-13.37 

2.10 

-21.00 

-8.25 

60 

3435 

-14.75 

-15.25 

2.18 

-23.00 

-7.50 

61 

1281 

-14.34 

-14.75 

1.95 

-22.50 

-9.00 

62 

8513 

-14.46 

-16.50 

6.28 

-36.75 

-2.00 

63 

11399 

-14.94 

-17.25 

5.78 

-37.50 

-2.75 

64 

10983 

-14.47 

-17.50 

5.40 

-34.75 

0.00 

65 

13363 

-13.51 

-15.50 

5.11 

-34.75 

-2.25 

66 

17312 

-11.74 

-12.75 

4.04 

-33.00 

1.75 

67 

3701 

-13.62 

-14.75 

2.99 

-27.25 

-5.00 
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Table  A4:  Backscattering  Coefficient  for  Mangrove  Forest  for  L-VV  POLSAR  Darwin 


N  Samples 

a°  (dB) 

Median  (dB) 

Std.  Dev.  (dB) 

Min  (dB) 

Max  (dB) 

27 

194 

-11.36 

-12.12 

2.87 

-19.50 

-6.00 

28 

99 

-10.76 

-11.00 

2.44 

-18.50 

-6.00 

29 

28 

-10.43 

-10.99 

2.55 

-15.75 

-5.50 

30 

142 

-11.73 

-12.75 

2.97 

-20.00 

-6.75 

31 

164 

-11.97 

-12.62 

2.39 

-18.75 

-8.00 

32 

52 

-12.50 

-12.75 

2.64 

-21.50 

-8.75 

33 

42 

-10.90 

-11.50 

2.86 

-19.75 

-5.25 

35 

46 

-10.69 

-11.50 

2.80 

-18.00 

-5.50 

40 

75 

-11.01 

-11.75 

2.30 

-17.75 

-5.00 

41 

29 

-12.08 

-13.00 

1.72 

-15.25 

-8.50 

42 

99 

-11.39 

-11.75 

1.98 

-17.75 

-8.00 

43 

105 

-10.34 

-11.25 

2.48 

-19.00 

-4.75 

44 

82 

-11.59 

-12.00 

2.58 

-18.25 

-6.25 

45 

114 

-10.61 

-11.12 

2.57 

-19.75 

-4.25 

46 

103 

-10.60 

-11.00 

2.31 

-18.25 

-6.50 

47 

178 

-10.04 

-11.00 

2.67 

-18.75 

-4.00 

48 

87 

-10.27 

-10.75 

2.20 

-16.50 

-6.50 

49 

47 

-9.80 

-10.50 

2.03 

-14.50 

-6.25 

50 

46 

-9.60 

-10.12 

2.15 

-14.25 

-6.25 

51 

44 

-9.24 

-10.25 

2.57 

-15.75 

-4.50 

52 

76 

-9.99 

-10.62 

1.93 

-15.00 

-6.50 

53 

55 

-10.34 

-10.75 

2.15 

-16.75 

-6.75 

54 

85 

-9.82 

-10.25 

2.39 

-18.75 

-3.25 

55 

87 

-8.40 

-9.00 

2.34 

-15.25 

-3.25 

56 

77 

-8.91 

-9.75 

2.79 

-16.50 

-2.75 

57 

41 

-8.19 

-9.25 

2.78 

-14.50 

-1.00 

58 

8 

-9.76 

-10.75 

2.41 

-14.50 

-7.00 

61 

28 

-8.65 

-9.86 

3.17 

-17.00 

-4.50 

62 

52 

-7.43 

-8.37 

2.64 

-17.00 

-2.00 

63 

51 

-7.70 

-8.25 

2.49 

-14.00 

-3.25 

64 

9 

-7.80 

-7.25 

2.68 

-13.50 

-5.25 

65 

41 

-7.90 

-9.25 

3.11 

-14.00 

-3.00 

66 

72 

-8.13 

-8.62 

2.30 

-13.25 

-3.00 

67 

50 

-8.93 

-9.25 

1.88 

-14.00 

-5.25 
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Appendix  B:  Semi-Empirical  Sea  Clutter 

Models 


B.l  Georgia  Institute  of  Technology  (GIT)  Sea  Clutter  Model 

[32,  1] 

(a)  Reflectivity  equations,  (dB  m2/m2) 

For  horizontal  polarisation, 

a°HH  =  101og[3.9  x  l0-6\ip0AGmGuGw]  (Bl) 

For  vertical  polarisation, 

a°vv  =  a°HH  -  1.73 In (ha  +  0.015)  4-  3.76  ln(A)  +  2.46 ln(^  +  0.0001)  4-  22.2  (B2) 

(b)  Interference  Factor 

a4 

Gm  =  (1  +  a4) 

(14.4A  4-  5 .5)ipha 
a  =  X 

Average  wave  height,  ha(m)  =  4.52  x  10  3V^'5 


(c)  Wind  Direction  Factor 

Gu  =  exp{0.2  cos  0(1  —  2.8V7) (A  +  0.015)  °'4} 


(d)  Wind  Speed  Factor 


Gw 

q 


'  1.9414, 

(14-1^/15.4) 

1.1 

(A  +  0.015)0-4 


where 

0  =  grazing  angle  (radians); 

A  =  radar  wavelength  (m); 
Vw  =  wind  speed  (m/ s) 
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B.2  Hybrid  (HYB)  Sea  Clutter  Model  [35,  1] 

(a)  Reflectivity  equations,  (dB  m2/m2), 

<t°  =  a°(ref)+Kg  +  Ks  +  Kp  +  Kd  (B3) 

(b)  Reference  Reflectivity  (dB  m2/m2)  for  ip  =  0.1°,  <p  =  0°,  VV  polarisation  and  sea 
state  5. 

£T° (re/)  =  24.4  log  /  -  65.2  (B4) 


(c)  Grazing  Angle  adjustments,  Kg  (dB) 

Reference  Grazing  Angle  (°),  ipr  =  0.1° 

Transitional  Angle  (°),  ipt  =  sin-1(0.0632A //irms)° 

RMS  Wave  Height  (m),  Hrms  =  0.031  .SS2  where  SS  =  Douglas  Sea  State 


(i)  For  ipt  >  ipr, 


for  Ip  <  1pr, 
for  Ipr  <  V’  < 

20 log  ^  +  10 log  Jj,  for  ipt  <ip  <  30°. 


'0, 

Kg=  <{  20  log 


(ii)  For  ipt  <  ipr, 


Kg  = 


( 0,  for  ip  <  ipr, 

\l0  log  for  ip  >  ipr- 


(d)  Sea  State  Adjustment  (dB),  Ks  =  5 (SS  —  5) 

(e)  Polarisation  Adjustment,  Kp  (dB), 

(i)  Average  Wave  Height  (m),  ha  =  0.08 SS2 

(ii)  Vertical  Polarisation,  Kp  —  0 

(iii)  Horizontal  Polarisation, 

Kp  =  1.71n(/ia  +  0.015)  —  3.8  In  A  —  2.5  In 

(f)  Wind  Direction  Adjustment,  K d  (dB) 

Kd  =  (2  +  1.7  log  (cos  <p  -  1) 


57.3 


+  0.0001  -  22.2,  (B5) 


(B6) 


where  <p  is  the  radar  look  angle  with  respect  to  the  wind  direction  (<p  —  0  for  upwind). 
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B.3  Technology  Service  Corporation  (TSC)  Sea  Clutter  Model 
[40] 

(a)  Low  Angle  Factor,  Ga 

Surface  height  standard  deviation,  az  =  0. 115. S'51'95  (B7) 

=  14‘l.  (B8) 

Ga  =  <£'/(! +  0  (B9) 

where 

xp  =  grazing  angle  (radians); 

A  =  radar  wavelength  (ft.); 

SS  =  Sea  state  (Douglas  sea  state); 
az  =  sea  surface  height  standard  deviation. 

(b)  Wind  Speed  Factor, 

Wind  velocity,  Vu 

Q 

Ax 

A2 

A3 

Ai 
A 


=  6.2.550'8  (BIO) 

=  <-6 


=  1  4-  0.35. <3 


2.63Ai 

A2A3A4 

'Vu  +  4.01 A 

15 
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For  vertical  polarisation  (dB), 

(jyy  =  a°HH  -  101og10  [1.73  ln(2.507cr2  +0.05)  +  3.76  In  A  +  2.46  ln(sin^  +  0.0001)  +  19.8] 

(B13) 


(d)  Specular  Returns 

For  horizontal  and  vertical  polarisation, 

(Jyy  —  & [JH  ~  +  COt  /30  exp 


tan2  9 
tan2  0o 


where 


9  =  n/2  —  ip 

10 log /i  =  -5  +  12.5(logA  -  log(0.5))dB 

f  10.1  +  1.65  55,  55  <2 

“  113.4  +  0.7(55-2),  55  >2 
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concept  of  a  two-scale  sea  surface  model  with  directional  sea  spectrum  and  simple  facet  specular  return 
at  near-normal  incidence.  Several  sea  clutter  models  have  been  collated  with  data  collected  by  the 
Naval  Research  Laboratory  four-frequency  radar  system  under  varying  sea  conditions,  and  polarimetric 
synthetic  aperture  radar  images  of  the  North-West  of  Australia.  The  results  of  the  comparison  have 
shown  that  the  composite  sea  clutter  model  better  predicts  the  backscattering  coefficient,  cr°,  of  sea 
surface  returns  for  varying  sea  conditions,  radar  look  directions,  and  grazing  angles.  In  addition,  a  ladar 
backscattering  model  for  bare  soil  showed  good  agreement  between  modelling  results  and  POLSAR  data 

for  returns  from  coastal  saline  flat  surfaces  surrounding  the  North-West  coastal  areas. _ 
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